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CLASSICAL  AND  GENERALIZED 


SOLUTIONS  OF  TIME-DEPENDENT  LINEAR 
DIFFERENTIAL  ALGEBRAIC  EQUATIONS^ 


BV 

Patrick  J.  Rabier  and  Werner  C.  Rheinboldt^ 

Abstract.  A  reduction  procedure  is  developed  for  linear  time-dependent  differential  alge¬ 
braic  equations  (DAEs)  that  transforms  their  solutions  into  solutions  of  smaller  systems  of 
ordinary  differential  equations  (ODEs).  The  procedure  applies  to  classical  as  well  as  distri¬ 
bution  solutions.  In  the  case  of  analytic  coefficients  the  hypotheses  required  for  the  reduction 
not  only  are  necessary  for  the  validity  of  the  existence  and  uniqueness  results,  but  they  even 
allow  for  the  presence  of  singularities.  Straightforward  extensions  including  undetermined 
systems  and  systems  with  non-analytic  coefficients  are  also  discussed. 


1.  Introduction. 

A  time-dependent  linear  differential-algebraic  equation  (DAE)  is  a  problem  of  the  form 
(1.1)  A{t)x  +  B(t)x  =  b{t),  t^J, 

where  J  C  R  is  an  open  interval,  G  £(R”)  and  A(t)  is  not  invertible  for  any 

t  ^  J.  In  the  simplest  case,  when  A(t)  is  diagonal,  this  condition  implies  that  -  at  least 
locally  -  some  of  the  scalar  equations  in  the  system  (1.1)  are  purely  algebraic,  which  justifies 
the  terminology"  “differential-algebraic  equation”.  In  order  to  provide  a  perspective  to  our 
work,  we  begin  with  a  brief  review  of  what  is  currently  known  about  such  problems  without, 
however,  including  some  earlier  contributions  now  subsumed  by  more  general  theories. 

A  significant  simplification  is  the  case  when  A{t)  and  B{t)  in  (1.1)  are  constant.  The 
corresponding  existence  and  uniqueness  theory  for  initial  value  problems  has  been  known 
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for  some  time:  A  fa!!  account  can  be  fou'u!  in  Griepentrog  and  Maerz  [GrMS6]  who 
credit  Gantmacher's  theory  of  matrix  pencils,  [G59].  The  book  [GrM86]  also  contains 
a  treatment  of  the  general  system  (1.1)  utilizing  a  condition  of  “transferability"  which, 
broadly  speaking,  cor-esponds  to  a  generalization  of  the  index  1  case  of  the  constant 
coefficients  theory.  Subsequent  generalizations  to  the  case  of  the  index  2  and  3  have  been 
obtained  by  Maerz  ([M89a],  [M89b]). 

A  different  approach  is  taken  in  Campbell  and  Petzold  [CP83]  who  consider  "solvable" 
systems  with  analytic  coefficients.  The  soh'ability  assumption  requires  that  an  existence 
and  uniqueness  theory  similar  to  that  for  linear  ODE's  is  available.  The  result  of  [CP83] 
is  that  such  systems  can  be  transformed  into  the  form 

(1.2a)  y  +  C(f)y  = /(#),  tej. 

(1.2i)  N{t)z  +  z  =  g{t),  iej. 

where  C(f)  and  J\'’(f)  are  analytic  functions  of  t  and  N{t)  is  nilpotent  upper  (or  lower) 
triangular  for  all  f  €  From  the  structure  of  N{t),  it  follows  that  the  operator  A'(f)^ 

n 

is  nilpotent,  so  that  (1.2b)  has  the  unique  solution  ^  ^  and  (1.2a)  is 

an  explicit  ODE.  But  no  calculable  criterion  for  solvability  is  given  in  that  paper. 

Later,  Campbell  [C87]  proved  results  of  similar  type  for  the  case  of  non-analytic  coef¬ 
ficients  in  (1.1)  and  provided  a  calculable  necessary  and  sufficient  criterion  for  solvability. 
This  condition  requires,  in  particular,  the  characterization  of  the  null-spaces  of  some  aug¬ 
mented  systems  for  each  t  ^  J .  He  also  showed  that  the  solutions  of  (1.1)  can  be  obtained 
as  the  solutions  of  an  explicit  ODE  in  R”  satisfying  consistent  initial  conditions. 

Recently,  Kunkel  and  Mehrmann  [KuMe92]  have  taken  a  different  approach.  Under 
various  constant -rank  conditions  they  show  that  the  solutions  of  (1.1)  can  be  obtained 
in  the  form  x  —  Uy  where  U  is  a.  smoothly  parametrized  family  of  linear  isomorphisms 
of  R"  and  y  =  (yi  -  y2,y3)'^  solves  a  simple  DAE  of  the  form  yi  =  ci(f),  yo  =  C2(t), 
0  =  C3(f),  where  c  =  (ci,  02,03)^  depends  upon  ^  and  its  derivatives.  Solvability  requires 
the  dimension  of  the  third  \'ariable  to  be  zero,  for  otherwise  y^  is  arbitrary  if  C3  =  0  or  no 
solution  exists  when  C3  7^  0,  and  hence  either  uniqueness  or  existence  fails  to  hold.  This 


approach  does  not  require  a  soh'ability  assumption,  for  the  dimension  of  the  third  \'ariable 
is  determined  in  the  process  of  finding  U .  On  the  other  hand,  when  the  DAE  is  solvable, 
the  hypotheses  made  in  [KuMe92]  to  reach  this  conclusion  are  sufficient  but  not  neccesary. 

FmaJly,  we  mention  new  theories  for  general  implicit  DAE’s  F{t.x,x)  =  0  developed  by 
the  authors  ([RRh91],  [RRh92]),  that  can  be  applied  to  hnear  problems.  But.  of  course, 
these  general  theories  do  not  capture  the  special  features  due  to  the  linear  structure. 

In  all  mentioned  results  only  classical  solutions  are  considered:  that  is.  all  solutions 
X  ■.  J  R"  of  (1.1)  are  assumed  to  be  (at  least)  differentiable.  As  far  as  we  know, 
the  problem  of  characterizing  genertilized  (distribution)  solutions  of  (1.1),  when  the  right- 
hand  side  6  is  a  distribution,  htis  been  successfully  investigated  only  for  constant  A  and 
B  and  in  special  cases  (see  [VLKaSl],  [Co82],  [Ge93]).  The  difficulty  with  distribution 
solutions  is  that  arguments  from  the  “classical”  theory  involving  the  numerical  \"alues  of 
the  solutions  can  no  longer  be  used.  This  is  the  reason  why  most  approaches  do  not  extend 
to  distributions.  It  appears  that  the  method  in  [KuMe92]  may  allow  for  such  an  extension, 
but  there  is  no  comment  to  that  extent  in  that  paper. 

In  this  paper,  we  focus  mainly  on  the  case  when  A  and  B  are  analytic.  Under  com¬ 
putationally  verifiable  assumptions  -  weaker  than  the  solvability  condition  of  [CP83]  -- 
we  prove  that  both  the  classical  and  distribution  solutions  of  (1.1)  are  of  the  form 
x(t)  =  r(f )ii,(f )  -f  f(  0  where  solves  a  smaller  system 

(1.3)  A^x^  +  B^x^  =  b^,  inj, 

and  .4^(f)  is  invertible  for  all  but  isolated  values  of  t.  It  turns  out  that,  when  (1.1)  is 
solvable  in  the  sense  of  [CP83],  (1.3)  can  also  be  derived  by  the  method  briefly  sketched 
at  the  end  of  that  paper.  But  extending  this  method  without  the  solvability  hypothesis 
requires  a  much  closer  investigation  of  its  mechanism,  especially  in  the  distribution  case. 

Equation  (1.3)  is  obtained  at  the  step  of  a  recursive  process  involving  a  number 
of  arbitrary  choices.  We  show,  however,  that  both  the  so-called  index  u  of  (1.1)  and  any 
singularities  of  are  independent  of  these  choices.  We  also  prove  that  invertibility  of 
A„(t)  for  each  t  ^  J  \s  necessary  and  sufficient  for  solvability  in  the  sense  of  [CP83].  But 
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since  solvability  is  not  a  prerequisite  for  the  validity  of  our  reduction,  singularities  may 
exist  in  (1.3)  and  the  methods  of  the  Fuchs- Frobenius  theorj"  are  available  for  the  study 
of  the  related  phenomena. 

We  also  outline  a  more  or  less  straightforward  generalization  of  our  reduction  procedure 
to  characterize  the  classical  solutions  of  any  analytic  system  (1.1).  although  either  existence 
or  uniqueness  (or  both)  for  corresponding  initicd  \^lue  problems  is  necessarily  affected 
whenever  the  hypotheses  of  the  main  case  break  down.  Distribution  solutions  can  be 
characterized  as  well  for  more  general,  but  perhaps  not  all  non-analytic  systems. 

The  inherent  globality  of  the  reduction  is  especially  important  to  hamdle  boundary  \'alue 
problems,  which  reduce  to  boundary  value  problems  for  (1.3). 

Finally  we  excunine  problems  (1.1)  in  the  case  when  A  and  B  are  sufficiently  smooth 
but  no  longer  analytic.  It  turns  out  that  then  additional  hypotlieses  are  needed  (except 
in  the  index  1  case)  to  compensate  for  the  lack  of  analyticity.  Because  of  this,  we  are  not 
able  to  cover  fully  the  “solvable  case”  of  [C87].  On  the  other  hand,  our  method  remzuns 
valid  with  weaher  smoothness  hypotheses,  say  C”,  for  ,4,  B  and  b  in  (1.1),  continues  to  be 
applicable  for  distribution  solutions,  zmd  can  be  easily  generalized  to  handle  undetermined 
problems  satisfying  a  Fredholm  alternative. 

This  work  and  [KuMe92]  are  certainly  close  in  spirit,  but  major  differences  exist  in  the 
execution  of  the  main  idea.  Our  reduction  appears  to  be  technically  much  simpler  and 
allows  for  clearer  proofs.  More  importantly,  it  is  this  relative  simplicity  that  permitted  a 
deeper  analysis  under  weaker  hypotheses. 

Our  approach  makes  crucial  use  of  Kato’s  “transformation  functions”,  and  the  next 
section  presents  a  brief  review  and  some  applications  of  this  concept.  Very  closely  related 
but  somewhat  less  precise  results  can  be  found  in  [D64]  and  [SiBTOj.  Throughout  the  paper 
analytic  mappings  will  be  referred  to  as  “mappings  of  class  C"”. 

2.  Transformation  functions. 

A  main  tool  for  the  proof  of  the  globality  of  the  reduction  procedure  is  the  concept 
of  “transformation  functions”  introduced  by  T.  Kato  [K50].  We  use  only  the  following 
simplified  version: 
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Theorem  2.1.  Let  J  C  R  be  an  open  interval  and  P  G  C*(J7:  £(R" )),  1  <  <  oc  or 

k  =  u),  a  family  of  orthogonal  projections  onto  VT(<)  =  rge  P{t).  t  G  J .  Then,  for  any 
fixed  to  £  J ,  there  exists  a  mapping  U  E  C'‘{J,0{n)),  where  0{n )  denotes  the  orthogonal 
group  of  R”.  such  that 

(2.1)  P{t)  =  U{t)P{to)U(tf,  Vf  ej. 

Hence  U(t)P{to)  is  a  linear  isomorphism  oflV(to)  onto  for  t  G  J . 

The  proof  of  a  more  general  result  involving  arbitrary'  projections  and  a  complex  \^riable 
in  place  of  t  may  be  found  in  [K82]  (p.  113).  It  will  be  useful  to  outline  the  proof  of  the 
somewhat  simpler  Theorem  2.1.  The  transformation  U  is  determined  as  the  solution  of 
the  linear  initial  value  problem 

(2.2)  U  =  [P,P]l\ 

(2.3)  Uito)  =  In, 

where  [P,  P]  =  PP  —  PP  is  the  commutator  of  P  and  P.  Because  of  P^  =  P.  the 
commutator  [P,  P]  is  skew  symmetric  whence  =  U^U  +  U^U  =  P^[P,P]P  - 

U'^[P,  P]U  =  0.  Therefore,  with  (2.3),  it  follows  that  U{t)^U{t)  =  and  hence  that 
U[t)  G  0{n)  for  i  G  J .  For  the  proof  of  (2.1)  we  use  P^  —  P  which  yields  PP  +  PP  =  P 
and  thus  PPP  =  0  and  [P,P]P  =  PP,  P[P,  P]  =  — PP.  With  this  we  find  for  F  =  PV 
that  Y  —  PU  +  PU  =  (P  —  PP)U  =  PPU  =  [P, P]y  which  shows  that  Y  solves  (2.2) 
with  the  initial  condition  Y{to)  =  P{to).  Since  the  general  solution  .Y  of  (2.2)  is  given  by 
A’  =  UX{to),  it  follows  that  Y  =  UPito)  and.  therefore,  that  PU  =  UP{to)  which  is  (2.1). 

For  later  use  we  observe  that  the  above  relation  Y  =  PPU  also  reads  Y  =  PY  and 
hence  that  Y  satisfies  the  differential  equation 

(2.4)  Y  =  PY. 

From  P^  =  P  we  infer  that  P^  =  P.  so  that,  by  (2.4),  Y^  satisfies  the  differential  equation 

(2.5)  Y^  =  r^P. 
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Our  applications  of  Theorem  2.1  rest  on  the  following  result  (see  also  Lemma  2.2). 

Lemma  2.1.  Let  J  C  R  be  an  open  interr'aJ.  n  >  0  an  arbitrary-  integer,  and  M  G 
C*(  £(R")),  1  <  A;  <  oo  or  A-  =  a;,  such  that  rank  M(i)  =  r.  Vt  Q  J.  Then,  the 

orthogonal  projections  onto  rge  M{t)  and  ker  M{t),  respectively,  are  cf  class  C*  in  J7. 

Proof.  Let  t  £  J  he  given  and  consider  an  orthonormal  basis  {ui,.. .  .Ur}  of  rge  M{t). 
Let  ,  Ur  be  linearly  independent  vectors  such  that  u,  =  M{t)v,,  1  <  i  <  r.  The 

vectors  Ui{t)  =  M{t)vi  remain  hnearly  independent  for  t  near  t.  say,  for  t  in  the  open 
interval  Ji  G  J  containing  t.  In  other  words,  since  rank  M{t)  =  r,  they  form  a  basis 
of  rge  M{t)  for  t  ^  and  they  are  of  class  in  since  this  is  true  for  M.  By- 
applying  the  Gramm-Schmidt  process  to  . .  ,Ur{t),  we  obtain  an  orthonormal  basis 

{tri(t),...  ,iVr{t)}  of  rge  M{t)  for  t  £  Jp  Since  the  Gramm-Schmidt  process  involves 
only  algebraic  operations,  the  vectors  lei(t),...  ,Wr{t)  are  ageiin  functions  of  t  E  Jp 
Because  the  orthogonal  projection  of  R”  onto  rge  M (t)  can  be  expressed  as  a  sum  of  dyadic 
products  of  wi{t), . . .  .  Wr{t),  it  follows  that  it  is  of  class  in  and,  hence  also  in  J 
since  t  €  J  was  arbitrary.  The  corresponding  result  for  the  orthogonal  projection  from  R" 
onto  k.exM{t)  follows  directly  by  replacing  M{t)  with  in  the  above  argument.  □ 

Theorem  2.2.  Let  J  C  R  be  an  open  interi'aJ.  n  >  0  an  arbitrary  integer,  and  M  G 
£(R")),  1  <  A-  <  oo.  such  that  rank  M{t)  =  r  for  all  t  £  J :  Then  the  following 
operators  exist: 

(i)  S  e  J;£(R^R"))  such  that  S(t)  €  GL(R^,  rge  M{t)),  Vt  E  J, 

(ii)  T  E  C'^(J;£(R"-^R"))  such  that  Tit)  E  GL(R"-^ ker  Af(t)),  Vt  E  J. 

(Hi)  V  E  C^(J;£(R",R''))  such  that  mo  ^  GLirge  M(t},R^),  Vt  E  J.  where 

V'  can  be  chosen  so  that  kerF(t)  =  [  rge  for  all  t  ^  J . 

Proof,  (i)  Denote  by  P{t)  the  orthogonal  projection  from  R"  onto  rge  A/(t),  t  E  J ■  By- 
Lemma  2.1,  we  have  P  E  C*(  J”;  £(R”)).  Consequently,  the  mapping  U  E  J":  £(R”)) 
of  Theorem  2.1  satisfies  U{t)P(to)  E  G£(rge  M(<o).  rge  A/(t)).  Vt  E  J.  Now.  choose 
5o  E  £(R''.R")  such  that  5o(R'')  =  rge  M{to)  and  set  Sit)  =  Uit)So.  Vt  E  J. 

(ii)  Let  n(t)  be  the  orthogonal  projection  from  R"  onto  kerA/(t).  t  ^  J .  By-  Lemma 
2.1,  we  have  II  E  Jf;  £(R" )).  Hence  the  corresponding  mapping  U  E  £(R”)) 
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of  Theorem  2.1  satisfies  L’’(On(to)  €  Gi(ker A/(to ), her j).  Vt  ^  J.  Now.  choose 
To  fc  £(R"~'',R")  such  that  To(R"“'')  =  kerA/(fo)  and  set  T{i)  =  U(t)To.  Vt  6  ^7. 

(iii)  With  P  and  U  as  in  the  proof  of  part  (i),  let  L  G  £(R",R'')  be  such  that 

L(rge  M(to))  =  R'’.  Then,  set  V'(t)  =  LP(to)U{t)^,  'it  6  J.  It  is  readily  seen  that 

V'(0|rge  Af(,)  ^  rge  Furthermore.  kerV'(0  =  ker  P{to)U{t)^  =  {x  e  R^  : 

U(t)'^x  e  [rge  A/(#o)]'^}  =  U{t)  [rge  A/(#o)]'‘'-  As  U{t)  e  0{n)  maps  rge  A/(to)  onto 

rge  A/(t),  it  also  maps  [rge  A/(^o)]'*‘  onto  [rge  A/ff)]-*-,  and  hence  kerV'(t)  =  [rge  A/(0]'^- 
it  6  J,  as  claimed.  □ 

Note  that  by  (2.4)  and  (2.5)  the  mapping  T  of  Theorem  2.2  was  found  to  be  the  unique 
solutions  of  the  initial  value  problem 

(2.6)  T  =  riT,  T(<o)  =  To. 

where  to  ^  <7^  is  some  fixed  value.  n(t)  denotes  again  the  orthogonal  projection  onto 
kerA/(t),  t  G  7,  and  To  €  £(R”“'',R")  is  any  operator  such  that  To(R"~'')  =  kerA/(to)- 
Analogously,  the  mapping  V  is  the  solution  of 

(2.7)  V  =  VP.  V{to)  =  LP{t,), 

where  P(t)  is  the  orthogonal  projection  onto  rge  M{t},  t  ^  J  and  L  G  £(R".R'')  is  any 
operator  such  that  T(rge  M{to))  =  R*". 

In  the  analytic  case,  the  assumption  that  rank  M{t)  =  r.  it  ^  J ,  can  be  dropped 
without  major  consequences  for  Theorem  2.2  because  of  the  following  mrlanL  of  Lciiiina 
2.1: 

Lemma  2.2.  Let  7  C  R  be  an  open  interval,  n  >  1  any  integer,  and  M  G  £(R")). 

Set  r  —  maxtgj  rank  M{t).  Then,  there  is  a  subset  S  C  J  of  isolated  points  such  that 
rank  M(t)  =  r  if  and  only  if  t  E  \  S.  Furthermore,  the  orthogonal  projections  onto 
rge  M(t)  and  kerM{t).  t  G  7  \  >S.  which  are  analytic  on  t  ^  J  \  S.  can  he  extended  as 
analytic  functions  over  the  entire  interval  J . 


Proof.  Since  rank  Mlt)  acliieves  only  the  discrete  values  0. 1.  ■  •  .n.  there  exists  a  f  ^  j 
such  that  rank  M(t)  =  r.  Identifying  M(t)  with  its  matrix  in  the  canonical  basis  of  K". 
we  see  that  some  r  x  r  minor  of  M{t)  does  not  vanish  identically  in  J .  and  hence,  by 
analyticity,  can  \'anish  only  on  some  subset  S  of  isolated  points  of  J .  Since  r^e  M(t)  = 
rge  ,  Vf  G  J,  we  may  replace  M  by  M without  modifying  rge  M{f)  or 

affecting  the  analyticity.  It  then  follows  from  the  symmetry  of  MM^  and  Theorem  G.l. 
p.l39  in  [KS2]  that  the  orthogonal  projection  P(t)  onto  rge  M(t).  t  ^  J  \  S.  can  be 
extended  as  an  analytic  function  in  J .  (In  the  notation  of  [KS2].  P  is  the  sum  of  the  eigen- 
projections  Ph  corresponding  to  the  eigenvalues  \h  o{  that  do  not  vanish  identically 

in  because  eigenprojections  of  symmetric  operators  are  orthogonal  projections).  Let 
now  P{t)  denote  the  orthogonal  projection  onto  kerA/(f).  t  ^  J’  \  S ■  For  the  proof  that 
P  can  be  extended  as  cin  analytic  fimction  in  J  it  suffices  to  replace  M(t)  by  Mit)^ M{t) 
and  to  use  once  again  Theorem  6.1.  p.l39  in  [K  82].  □ 

Let  M,  r,  and  S  be  as  in  Lemma  2.2.  and  call  P(t)  the  orthogonal  projection  onto 
rge  M{t),  t  G  \  S.  so  that  P  G  £(R”)).  By  the  denseness  oi  J  \  S  \r\  J  and  the 

relation  PM  =  M  in  J"  \  5,  we  find  that  PM  =  M  in  J  and.  in  particular,  rge  M(t)  C 
rge  P{i)  for  t  G  S.  Moreover.  P{t)  remains  an  orthogonal  projection  at  points  of  5  as 
follow  from  the  relations  P^  =  P  and  P^  =  P  extending  from  J  \S  to  J .  In  other  words. 
P{t)  is  a  projection  onto  a  subspace  containing  rge  M{t).  Let  now  Pj  G  £(1R" ) )  be 

any  other  family  of  projections  such  that  rge  Pi(t)  =  rge  M{t),  Vt  G  S  \  S .  Obviously,  in 
J  \  S,  and  hence  in  J.  we  have  PPi  =  Pi  so  that  rge  P\{t)  C  rge  P{t).  Vf  G  J ■  By  an 
exchange  of  the  roles  of  P  and  Pi ,  it  follows  that  rge  Pi{t)  =  rge  P(/).  Vt  G  J .  This  shows 
that  even  at  points  of  S.  the  range  of  P{t)  is  independent  of  the  specific  choice  of  P.  For 
this  reason,  we  call  rge  P{t)  the  extended  range  of  M(t)  to  be  denoted  by  ext  rge  M{f). 
Then,  we  have 

(2.8a)  rge  M{t)  C  ext  rge  M(t),  Vt  G  J ~ 

(2.86)  rge  M(t)  =  ext  rge  M(t),  Vt  G  J  \  S . 

Let  n(t)  again  denote  the  orthogonal  projection  onto  ker.\/(t).  t  G  J  \  S.  so  that 
n  G  C'^{J\  £(IR")).  The  denseness  ol  J  \  S  in  J  and  the  relation  .Mil  =  0  in  \  5  imply 
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that  ,\/n  =  0  in  J .  In  particular,  for  t  €  5.  wp  have  rge  IKt)  C  ker.\/(f)  and  hence  lilt) 
projects  onto  a  subspace  contained  'n  kci  M{  t).  Let  now  Hi  t  C"'(J:  £(  IR"  ))  be  any  other 
family  of  projections  such  tha.  ge  Hil/)  =  ker.\/(n.  VY  €  \  5.  (3bviously.  mii  =  11 1 

and  riin  =  n  in  '  5.  and  hence  in  J.  so  that  rge  Odf)  =  rge  fljt).  vY  *=  J .  Thus,  the 
range  of  flff)  remains  independent  of  the  specific  choice  of  IT  at  points  of  S.  and  we  call 
rge  n(f)  th"'  restricted  null-space,  of  M{f)  to  be  denoted  by  rest  ker  M{t).  Then 

(2.9a)  kerd/(f  )  D  rest  ker  M{t).  Vf  € 

(2.96)  ker.\/(f)  =  rest  ker  .1/(0,  Vt  G  J/ \  5. 

By  Theorem  2.1.  ext  rge  M{t)  and  rest  ker  M(t)  arc  isomorphic  to  ext  rge  M{to)  and 
rest  ker  .1/(0).  respectively,  for  any  t.to  e  J.  Hence,  by  (2.Sb)  anrl  (2.9b).  we  have 
dim  ext  rge  M(t)  =  r.^t  e  J  and  dim  rest  ker  .1/(0  =  n  -  r.  Vf  G  J  which  shows  that 
the  inclusions  (2.Sa)  and  (2.9a)  are  strict  for  f  €  «S. 

Using  the  concepts  of  an  extended  range  and  restricted  null-space,  and  by  a  proof  similar 
to  that  of  Theorem  2.2.  we  obtain  now  the  following  version  of  that  result: 

Theorem  2.3.  Let  J  C  R  be  an  open  interxnJ.  n  >  0  any  integer,  and  .1/  G  ;  £(R"  ) ). 

Set  r  =  maxfgy  rank  .1/(0  sne?  let  S  C  be  the  set  of  isolated  points  where  rank  .!/( t)  <  r 
(see  Lemma  2.2).  Then  the  following  operators  exist: 

(i)  S  G  such  that  S{t)  G  GI(R’'.  e.\-f  rge  .!/(/) ).  Vt  G  J.  and.  in 

particular.  5(0  ^  GZ(R''.rge  M{t)).  '^t  ^  J  \  S. 

(li)  T  G  C-^(J;£(R"-MR"))  such  that  T{t]  G  GI(R"-0rest  her  .1/(0).  V/  G  J.  and. 
in  particular.  T(t)  €  GL(R"~''.  ker  .1/(0).  Vt  G  \  5. 

(Hi)  U  G  C'*'( J7;  £(R". R'"))  such  that  V(t)  G  GZ(rge  .!/(/), R'").  Vt  G  J.  and.  in  par¬ 
ticular.  U(0  €  GL{rs:e  M{t).R^).  €.  Lf  \  S  and  1'  can  be  chosen  so  that  kerU(0  = 

[ext  rge  .l/(f  )]-*-,  \/f  ^  J . 

Remark  2.1.  For  the  validity  of  (2.6)  and  (2.7),  n(#)  and  P{f)  must  be  chosen  as  the 
orthogonal  projections  onto  rest  ker  .1/(0  ''^nd  ext  rge  .1/(0.  respectively.  □ 
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3.  The  reduction  of  analytic  pairs. 

Let  J’  C  again  be  an  open  interval  and  A.  D  E  C'^(J ;  £(R” )).  We  denote  by  A~-  B  € 
C"'(.7;£l'R"  X  the  parametrized  family  of  operators 

(3.1)  (p.x)€R"xR"  ^  + V/eJ. 

Definition  3.1.  The  pair  (.4.  B)  is  called  regular  in  J  if  the  rank  condition 

(3.2)  rank  .4(f)  ►r- B(t)  =  n,  'it  ^  J . 

holds. 

For  later  use,  we  note  that  (3.2)  is  equivalent  with  either  one  of  the  conditions 

(3.2')  rank  [A{t)A(tf  +  B{t)B{tf]  ==  n,  VteJ. 

(3,2")  ker.4(f)^nker5(f)^  =  {0},  it  e  J. 

(3.2'")  rge.4(f)  +  rgeB(0  =  K",  itej. 

although  (3.2'")  will  not  be  particularly  useful  here. 

Suppose  that  (.4.  B)  is  regular  in  the  sense  of  Definition  3.1.  By  Lemma  2.2.  it  follows 
that 

(3.3)  rank  .4(f)  =  r, 

where  S  C  J  consists  of  isolated  points,  and  that  there  is  a  family  of  projections  P  G 
C"’{J\  C{W^))  with  rge  P{t)  =  ext  rge  .4(f),  Vf  €  J-  Set  Q  =  I  —  P.so  that  Q(t)  projects 
onto  a  complement  of  rge  .4(f),  Vf  6  J\S.  Clearly,  rge  .4(f)0J3(f)  =  rge  A{t)BQ(t)B{t)  = 
rge  ,4(f)0rge  Q(t)B{t).  Hence,  it  follows  from  (3.2)  and  (3.3)  that  rank  Q{t)B{t)  =  n  —  r. 
'it  ^  J  \  S.  or.  equivalently,  that 

(3.4)  dimker  Q(f)S(f)  =  r,  it^J\S. 
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Remark  3.1.  Using  (3.2)  we  see  easily  that  (3.4)  actually  holds  for  t  ^  J.  Indeed, 
given  u'  G  K",  there  axe  u,v  G  R"  such  that  i(t)u  4-  B[t)v  =  w.  Since  QA  —  0  in  J . 
it  follows  that  Q{t)B(t)v  =  Q{t)w.  Letting  u'  run  over  the  (n  —  r)-dimensional  subspace 
rge  Q{t)  we  infer  that  rank  Q(t)B(t)  >  n  —  r,  \/t  £  J .  Thus,  equality  holds  since,  by  (3.4) 
and  the  lower  semi-continuity  of  the  rank,  rank  Q{t)B(t)  <  n  —  r .  'it  ^  J .  As  a  result, 
"rest  ker  Q{t)B(ty'  could  be  replaced  by  “keTQ(t)B{ty'  in  the  subsequent  considerations. 
However,  we  have  chosen  to  ignore  this  fact,  as  this  will  later  enable  us  to  ascertain  that 
the  results  of  this  section  can  be  generalized  with  the  same  proofs,  to  a  case  when  (3.2) 
fails  to  hold.  □ 

Now.  let  C  G  C"'(J';£(RMR"))  be  such  that 

(3.5)  C(0  e  GI(RU rest  ker  (5(f)R(<)),  itej. 

The  existence  of  C  is  ensured  by  (3.4)  and  Theorem  2.3  (ii)  with  M  =  QB.  Analogously, 
let  D  G  C‘^(J;£(R".R'‘))  be  such  that 

(3.6)  D(t)  G  GI(ext  rge  A(t),R’'),  Vf  G  J, 

where  now  the  existence  of  D  is  ensured  by  (3.3)  and  Theorem  2.3  (iii)  with  M  =  .4. 
Remark  3.2.  For  t  E  J  \  S  we  have  keTQ{t)B{t)  =  {x  G  R"  :  B(t)x  G  rge  .4ft)}. 
which  shows  that  keTQ{t)B(t)  is  independent  of  the  particular  choice  of  Q;  that  is.  of  P. 
Consequently,  C(t)  is  independent  of  the  specific  choice  of  P  for  t  E  J  \  S  and  hence  for 
<  G  J.  □ 

Definition  3.2.  For  the  regular  pair  (A.B)  in  J  let  C  ^  C"'( JT";  £(R''.R" ))  and  D  G 
C'^iJ :  £(R”,  R'"))  satisfy  ( 3.5)  and  (3.6).  respectively.  Then  the  pair  (.4] .  B\ )  defined  by 

(3.7)  .4i,Fi  G  C‘^(J;£(R’')),  Aj  =  DAC.  Bj=D(BCaAC). 
is  a  reduction  of  (.4.  B)  in  J . 

Clearly,  a  reduction  (.4i,i5i)  of  a  regular  pair  {A.B)  is  not  unique,  since  it  depends 
upon  the  choice  of  C  eind  D.  However,  any  two  reductions  of  (A.B)  are  equivalent  under 
the  following  equivalence  relation  already  used  in  [CS7]  and  rediscovered  in  [KuMe92]. 


11 


Definition  3.3.  For  any  .4,  B,A,B£  £(R" ))  the  pairs  [A.  B)  and  (A,  B)  are  called 

equivalent  in  J,  to  be  denoted  by  (A,B)  ~  (A,B)  in  J ,  if 

(3.S)  A  =  MAN  and  B  =  M{BN  +  AN), 

for  some  M,  N  e  GI(R")). 

Lemma  3.1.  The  relation  of  Definition  3.3  is  an  equivalence  relation. 

Proof.  In  this  proof,  the  interval  J  will  be  implicit;  that  is,  {A,B)  ~  {A,B)  shall  mean 
that  this  relation  holds  in  JT.  Clearly,  (-4,5)  ~  {A,B)  is  trivial.  If  (-4.5)  ~  (-4.5) 
holds  then  -4  =  M~'^ AN~^  and  we  obtain  5  =  M~^{BN~^  +  -4^A^“^).  Indeed,  we  have 
M-^iBN-'^+Aj-^N-^)  =  {BN  +  AN)N-^  +ANj-^N-^  =  B  +  AiNN-^+Nj^N-^)  =  5. 
since  NN~^  +  N-^N~^  =1  =  0.  This  shows  that  {A,  5)  ~  (.4, 5). 

Now,  suppose  that  (^1,5)  ~  (-4,5)  and  (A,  5)  ~  (-4,5);  that  is,  A  =  MAN,  5  = 
M{BN  +  AN)  and  .4  =  MAN,  5  =  A/(5iV  +  AN)  for  M,N,M,N  €  GI(  J;  GI(R")). 
Then,  by  a  straightforward  calculation,  we  find  A  =  MMANN,  5  =  MN{BNN  + 
A-^{NN))  and  therefore  (.4,5)  ~  (-4,5).  □ 

In  the  proof  of  part  (iii)  of  our  next  theorem  and  in  other  places  later  on  we  require  the 
following  lemma. 

Lemma  3.2.  Let  P  ^  C‘^(J7^;£(R"))  be  a  family  of  projections  such  that  rge  P(t)  = 
ext  rgeA(t),  \ft  £  J  and  set  Q  =  I  —  P.  Then  for  any  integer  m  >  0  and  K  6 
C°(j7’;  £(]R'", R"))  we  have 

(i)  rge  K{t)  C  rest  ker  Q{t)B{t),  Wt  €  J , 
if  and  only  if 

(ii)  rge  B{t)I\{t)  C  ext  rgeA{t),  Vt  G  J , 

Note:  The  eqmvalence  between  (i)  euid  (ii)  above  is  not  true  for  fixed  t  £  J ,  although  (i) 
implies  (ii)  in  this  case  as  well. 

Proof.  Suppose  that  (i)  holds,  so  that  rge  K{t)  C  kerQ(<)5(t),  Vt  G  J .  This  is  equivalent 
with  QBK  =  0;  that  is,  rge  B{t)K{t)  C  kerQ(t)  =  rge  P{t)  =  ext  rge-4(<),  Vt  €  J,  which 
is  (ii).  Conversely,  assume  (ii)  to  hold,  so  that  QBK  =  0  and  thus  rge  K{t)  C  ker  Q{t)B{t). 


12 


6  J-  Let  n  €  C’^(j7;  £(R"))  be  the  orthogonal  projection  onto  rest  ker  Q{t)B{t). 
Since  rest  ker  Q{t)B(t)  =  ker  Q{t)B(t)  except  at  isolated  points  of  J,  we  find  that  (I  — 
n(f))/v(f)  =  0  except  at  isolated  points.  But  by  denseness  and  continuity,  it  follows 
that  necessarily  (I  —  IljA'  =  0;  that  is.  K  =  IIA',  whence  rge  K{t)  C  rge  IKf)  = 
rest  ker  Q{t)B{t),  'it  ^  J  which  is  (i).  □ 

Theorem  3.1.  Suppose  that  (A,B)  ~  (A, B)  in  J .  Then. 

(i)  rank  A{t)  =  rank  A{t),  'it  €  J. 

(ii)  (A.  B)  is  regular  in  J  if  and  only  if  {A,  B)  is  regular  in  J . 

(Hi)  If  (A,B)  and  (A,  B)  are  regular  in  J ,  any  reduction  (A],Bi)  of  (A.B)  in  J  is 
equivalent  to  any  reduction  (Ai,Bi)  of  (A.B)  in  J . 

Proof.  By  (3.8),  (i)  is  trivizJ.  To  show  that  (3.2),  or  equi\’alently  (3.2"),  either  fails  or 
holds  simultaneously  with  A  and  A,  it  suffices  to  show  that  if  (3.2")  is  valid  for  the 
pair  (.4,B),  then  it  also  holds  for  the  pair  (A.B).  Indeed,  symmetry’  of  the  equivalence 
relation  implies  that  the  converse  is  true  as  well.  Suppose  then  that  (3.2")  is  vedid  and  let 
t  ^  J  and  X  6  kerA(f)^  nkerB(f)^.  By  (3.8),  this  means  that  N{t)'^ A(t)'^M{t)^x  =  0 
and  N{t)'^B(t)^M{t)^x  +  N{t)^A(t)'^M{t)'^x  =  0.  Because  of  N'^(t)  G  (jX(R").  the 
first  relation  provides  that  M{t)^x  €  kerA(f)^.  and  hence  the  second  one  reduces  to 
N{t)^B{t)^M{t)'^x  =  0;  that  is,  M(f)^x  G  kerB(f)^.  Thus.  M{t)'^x  G  kerA(f)^  H 
kerB(f)^,  whence  M{t)'^x  =  0  by  (3.2").  Because  of  M{t)^  G  GAfR")  we  find  that  x  =  0 
which  proves  (ii). 

To  prove  (iii),  we  wite 

(3.9)  Ai  =  DAC,  Bi  =  DiBC  +  AC),  Aj  =  DAC.  Bj  =  D(BC  +  AC) 

where  C,  C  G  C‘^( JT":  £(R'',  R"))  and  D.D  G  C‘^(j7;£(R”,R’’))  satisfy  the  conditions  of 
Definition  3.2.  Note  here  that  the  integer  r  is  the  same  for  both  pairs  (A,B)  and  (A.B) 
by  part  (i)  of  this  theorem.  From  (3.8),  we  infer  that 

(3.10)  Ai  =  DMANC.  B,  =  DM{BNC  +  A^(.VC)). 
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Now  the  proof  hinges  on  the  remark  that  there  are  mappings  H.  L  ^  GL(W))  such 

that 

(3.11)  NC  =  CL 

(3-12)  -^(OAfiext  rge  A(<)  =  -^(O-^lext  rge  ^(0'  V/  €  JT", 

where  the  existence  of  H  and  L  will  be  established  further  below.  For  the  time  being,  we 
show’  how  the  validity  of  (3.10),  (3.11)  provides  for  the  proof  of  (iii).  From  (3.10),  (3.11) 
and  (3.12),  we  obtain  at  once 

(3.13)  ii  =  HDACL  =  HAiL, 

emd  Bi  =  DMBCL  +  HDA-^{CL).  Recall  that  C{t)  maps  into  rest  ker  Q{t)B{t)  for 
t  €  J7,  where  Q{t)  is  a  projection  onto  a  complement  of  ext  rge.4(t).  Therefore,  by  Lemma 
3.2  with  K  =  C,  B{t)C{t)  maps  into  ext  TgeA{t),  t  ^  J,  and  hence  (3.12)  implies  that 
DMBC  =  HDBC.  As  a  result,  B^  =  HDBCL  +  HDAf^{CL)  =  HD{BC  +  AC)L  + 
HDACL.  Using  (3.9)  we  thus  get  Bi  =  H{BiL  +  AiL)  and,  together  with  (3.13),  this 
implies  the  equivalence  of  {Ai,Bi)  and  {Ai,Bi)  since  H,L  £  C{W)). 

To  complete  the  proof,  we  now  verify  the  vedidity  of  (3.11)  and  (3.12).  By  hypothesis. 
B  =  M{BN  +  AN)  cmd,  by  construction  of  C,  B{t)C{t)  maps  into  ext  rge  A{t)  for  t  £  J . 
Since  ext  rge  A{t)  =  ext  rge  M{t)A{t)N{t)  =  ext  rge  M{t)A{t)  =  M(t)(ext  rge  .4(t)),  this 
amounts  to  saying  that  M {t){B {t)N {t)+ A{t)N{t))C{t)  maps  into  M{t){ext  rge  -4(t));  that 
is,  B{t)N{t)C{t)+A(t)N{t)C{t)  maps  into  ext  rge  A{t).  But  then,  B(t)N{t)C{t)  maps  into 
ext  rge  A(t).  By  Lemma  3.2,  N{t)C(t)  maps  into  rest  ker  Q{t)B{t)  =  rge  C{t).  'it  £  J 
with  Q  £  [J  \  C{'^'^))  as  in  (3.4)  and  (3.5).  By  reink  C{t)  =  r  and  N{t)  £  GX(R") 
we  have  rge  N{t)C{t)  =  rge  C{t).  On  the  other  hemd,  injectivity  of  C{t)  implies  that 
C{t)'^C{t)  £  GL{W),  it  £  J.  Set  Lit)  =  iC{t)Cit)'^)-^C{t)'^ Nit)C{t),  for  t  £  J, 
so  that  L  £  C‘^(j7:  £(R'^)).  In  fact,  L  £  C'^{J\GLiW)),  for  as  was  just  seen  above. 
N{t)C{t)  maps  onto  rge  C{t)  and  C{t)^  is  one-to-one  in  rge  C{t).  Thus,  rank  L{t)  =  r 
as  desired.  Next,  for  ^  €  R'"  we  have  L{t)^  =  (C(t)^C(t))~'(C(t)^JV(t)C(f)^,  whence 
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C{t)^C{t)L{t)^  =  C(t)^ N{t)C{t)^.  Since  C{t)^  is  one-to-one  in  rge  C(t)  =  rge 

this  gives  C{t)L(t)^  =  N{t)C{t)^-,  that  is,  CL  =  -VC  which  completes  the  proof  of  (3.11). 

To  prove  (3.12)  let  S  €  C'^( £(1R’',R'* ))  be  such  that  S(t)  6  GT(R''.ext  rge  A{t)). 
'it  €  J,  as  provided  by  Theorem  2.3  (i).  Note  that  D{t)S{t)  6  GZ(R'').  it  6  J7,  and 
set  H{t)  =  D{t)M{t]S{i){D{t)S{t)}~^ ,  so  that  H  €  C‘^(  J";  £(R’‘)).  In  fact,  we  obtain 
H  €  C^{J';  GL{R^))  because  M{t)  G  G£(R")  implies  that  M(t)S{t)  is  a  bijective  mapping 
into  M(i)[ext  rge  -^(t)]  =  ext  rge  A{t),  whereas  D{t)  is  one-to-one  in  ext  rge  A(t).  Now. 
for  X  G  ext  rge  A{t).  we  have  x  =  S(t)^  for  some  ^  G  R'^,  whence  D{t)x  =  D{t)S{t)^  and 
H{t)D{t)x  =  H{t)D{t)S{t)i.  Thus,  H{t)D{t)x  =  D{t)M(t)S{t)i  =  D(t)M{t)x.  which 
proves  (3.14).  □ 

Theorem  3.1  makes  it  easy  to  define  inductively  a  reduction  procedure  which,  up  to 
equivalence,  is  independent  of  any  particular  choice  that  must  be  made  at  each  step. 
Specifically,  suppose  that  the  pair  (A.B)  is  regular  in  JT^,  and  let  (Ai,j5i)  and  (,4i,.Bi) 
be  any  two  reductions  of  (A,B)  in  J.  By  Theorem  3.1,  {Ai,Bi)  ~  {A\,Bi)  cmd  hence 
( .4i ,  jBi  )  and  (A\,Bi)  are  simultaneously  regular  or  not  regular  in  JT".  If  they  are.  then  once 
again  by  Theorem  3.1.  any  reduction  {A2,B2)  of  (^i,  j3i)  is  equivalent  with  any  reduction 
(-42,52)  of  (-4i,.Bi).  Therefore,  (-42,52)  and  (^42, 52)  will  be  simultaneously  regular,  and 
so  on.  This  suggests  the  following  definition: 

Definition  3.4.  The  pair  (-4,5)  is  said  to  be  completely  regular  in  J  if  the  above  in¬ 
ductive  reduction  procedure  can  be  continued  indefinitely.  In  particular,  the  complete 
regularity  of  a  pair  [A.  B)  in  J  implies  its  regularity  in  J . 

-4  trivial  but  important  case  is  when  (-4,5)  is  regular  and  max^gj-  rank  A(t)  =  n  and 
thus  r  —  n.  Then  we  have  ext  rge  A(t)  =  R",  it  G  J.  whence  P{t)  =  /  is  the  only 
possible  choice  for  a  parametrized  family  of  projections  onto  ext  rge  A{t).  As  a  result. 
Q  =  I  —  P  =  0  and  therefore  QB  =  0  and  rest  ker  Q{t)B{t)  =  R",  it  G  J.  This 
shows  that  C(t)  =  D{t)  =  /  is  a  possible  choice  for  the  reduction  of  (-4,5)  for  which 
then,  by  (3.7),  (-4i,5i)  =  (-4,5).  Hence  it  follows  that  (-4,5)  is  completely  regular  and 
that  all  the  possible  reductions  are  equivalent  to  (.4,  5).  In  other  words,  up  to  equivalence, 
reduction  has  no  effect  on  the  pair  (-4, 5).  The  importance  of  this  case  is  that  it  will  always 
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occur  eventu2dly  at  some  stage  of  the  reduction  procedure  of  any  completely  regular  pair. 
Therefore,  despite  the  fact  that  the  reduction  can  be  theortically  continued  indefinitely,  it 
ceases  to  have  any  effect  (up  to  equivalence)  after  a  finite  number  of  steps.  This  is  the 
content  of  the  following  result: 

Theorem  3.2.  For  a  completely  regular  pair  {A,B)  in  J  set  >lo  =  -"i'  Bq  =  B  and  consider 
a  sequence  (Aj,  Bj),  j  =  0, 1, . . . ,  where  (AjBj)  is  some  reduction  of  {Aj-i,  Bj-\ )  in  J 
for  each  j  >  0.  Then. 

(3.14)  Tj  =  maxi^j  rank  Aj(t),  Vj  >  0, 

is  independent  of  the  specihc  choice  of  (Aj,Bj),  j  >  0,  and  with  r_i  =  n  we  have 

(3.15)  Aj{t)  €  £(R’'^-‘ ),  Vt  €  J,  V;  >  0, 
and 

(3.16)  r_i  =n  >  To  >  •••  >rj  >  •••  (>  0). 

Hence  there  is  a  smallest  integer  0  <  v  <  n  such  that  and  v  has  the  property 

(3.17)  .4^(<)  €  GF(R^-'-‘)  (=GF(R’'‘')), 

where  S„  G  J  consists  only  of  isolated  points,  and  v  is  also  characterized  by  the  fact  that 
it  is  the  smallest  integer  j  for  which  Aj{t)  G  (tF(R''''  ”'  )  for  some  t  ^  J . 

The  proof  is  a  trivial  consequence  of  Theorem  3.1  and  the  various  definitions. 

From  the  remarks  preceeding  the  theorem,  it  follows  by  (3.17)  that,  up  to  equivalence, 
the  reduction  has  no  effect  on  the  pair  (A^,B^)  and  hence  the  process  could  as  well  be 
stopped  at  this  point. 

Definition  3.5.  For  a  completely  regular  pair  (A,B)  in  J  the  integer  0  <  u  <  n  of 
Theorem  3.2  is  called  the  index  of  the  pair  in  J .  Since,  v  is  independent  of  th  '  choice  of 
the  reduction  at  each  step  of  the  process  the  in^'f^x  depends  only  on  the  pair  (A.  B). 
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4.  Global  reduction  of  differential-algebraic  equations:  classical  solutions. 

We  apply  now  the  reduction  procedure  of  the  previous  section  to  the  differential  algebraic 
equation 

(4.1)  A(t)i B{t)x  =  bit),  tej,  A,B  eC^iJ-XiR^)) 

cind  show  that  it  reduces  (4.1)  to  cui  explicit  linear  ODE  on  the  interval  JT.  Naturally  .4 
cuid  B  in  (4.1)  shall  now  be  called  the  coefficients  of  the  DAE. 

Definition  4.1.  The  DAE  (4.1)  with  anedytic  coefficients  is  reducible  (resp.  completely 
reducible)  if  the  pair  (A,B)  is  regular  (resp.  completely  regular)  in  J.  For  completely 
reducible  (4.1)  the  index  of  the  pair  (A,B)  in  J  is  the  index  of  (4.1 ). 

Theorem  4.1.  For  the  reducible  DAE  (4.1)  with  analytic  coefficients  let  (Aj ,  )  be  any 

reduction  of  the  pair  (A,  B)  in  J ;  that  is, 

(4.2)  Ai=DAC,  Bi=D(BC  +  AC), 

where  C  €  C‘^(j7;  £(1X'',  R")),  D  €  C‘^(J7’;£(R”,K’‘))  satisfy  (3.5)  and  (3.6),  respectively, 
and  r  =  maxt^j  rank  A(t).  Suppose  that  b  €  C*(47;R")  where  I  <  k  <  oc  or  k  =  tj,'. 
Then 

(i)  There  exist  uq  t  C*^(j7’;R")  such  that 

(4.3)  B(t)uo(t)  -  b(t)  e  ext  rge  A(t),  'it  £  J , 
and  indeed  (4.3)  holds  for 

(4.4)  uo  =  B'^iAA'^  +  BB^)-'b. 

(ii )  A  differentiable  mapping  x  -.  J  — >  R"  solves  (4.1)  if  and  only  if  for  given  uq  satisfying 
(4.3)  we  have 

(4.5)  X  =  Cxi  +  uo 
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where  xi  ■.  J  — ►  R’"  is  a  differentiable  solution  of  the  DAE 


(4.6) 


Ai(t)xi  -h  Bi(t)x  =  bt(t),  i  e  J, 


with  bi  given  by 


(4.7) 


bi  =  D(b  —  Buo  —  Atio). 


Proof,  (i)  By  the  reducibility  of  {A,B)  and  (3.2'),  we  have  A{t)A{t)^  +  B{t)B{t)^  ^ 
GZ(R”),  Vt  €  J.  Thus,  with  uq  given  by  (4.4)  and  uq  =  A'^{AA^  +  BB^)~^b,  we  obtain 
Avq  +  Buq  =  b  which  shows  that  (4.3)  holds. 

(ii)  Let  X  •.  J  — >  R”  be  a  differentiable  solution  of  (4.1)  so  that  B[i)x{i)  —  b{t)  £ 
rge  A{t),  V<  €  J.  For  uq  £  C*'(J7’;R")  chosen  according  to  part  (i)  we  have  B(t)(x(t)  — 
uo(f))  €  ext  rge  A(t),  Vf  £  J .  Then,  by  applying  Lemma  3.2  with  Q{t)  as  stated  there 
cind  K  £  C°(J';£(R^R"))  defined  by  K{t)  =  x{t)  —  ■uo(f)>  we  obtain  x{t)  -  uo(t)  £ 
rest  ker  Q{t)B{t).  But  now  (3.5)  implies  that  x{t)  —  uo(t)  =  C(t)xi(t)  for  a  (unique) 
ii(f)  £  R*",  and  xi(f)  =  [C(f)^C(f)]“^C(f)^(x(<)  —  uo(f))  shows  that  xi  :  >  R*" 

is  differentiable.  By  differentiating  x  =  Cxi  +  uq  we  obtain  ACxi  +  (BC  +  -4C)xi  = 
b  —  .4uo  —  Buq  and  hence,  after  multiplication  by  D,  that  Xi  solves  the  DAE  (4.6). 

Conversely,  let  X]  :  ^  R*"  be  a  differentiable  solution  of  (4.6).  Then,  by  (4.2),  we 

have 

(4.S)  0  =  D[-4(Gxi  +  Gxj  +  lio)  4-  B{Cx\  +  wq)  —  6]. 

By  (3.5)  and  Lemma  3.2  B{t)C{t)  maps  into  ext  rge  A{t)  while,  by  (4.3),  B{t)uQ{t)~b{t)  £ 
ext  rge  .4(f).  This  shows  that  the  bracketed  term  in  (4.8)  belongs  to  ext  rge  A(f)  for  each 
t  £  J .  But,  by  (3.6)  D{t)  is  an  isomorphism  of  ext  rge  A{t)  onto  R'",  Vf  £  J .  and  hence 
this  bracketed  term  vainishes.  Therefore,  x  defined  by  x  =  Cxi  -I-  uq  is  a  solution  of 
(4.1).  □ 


18 


Corollary  4.1.  Let  the  DAE  (4.1)  with  analytic  coefficients  be  completely  reducible  on 
J  with  index  i/  >  0  and  assume  that  b  £  with  v  <  k  <  oc  or  k  =  Set 

.4o  =  A.  Bo  —  B  and  consider  a  sequence  (Aj.Bj).  j  =  0.1 . where  (AjBj)  is  some 

reduction  of  (Aj-i.  Bj-i)  in  J  for  each  j  >  0;  that  is, 

(4.9)  Aj  =  Dj^iAj-iCj-i,  Bj  =  Dj-i(Bj^iCj-i  Aj-iCj-i) 

where  Cj-i  and  satisfy  the  condition  of  the  reduction  procedure  and  r^,  j  >  0  are 

defined  by  (3.16)  and  r_i  =  n.  Then  the  following  holds: 

(i)  There  exist  sequences  Uj,bj  £  ),  0  <  j  <  such  that  bo  =  b  and 

(4.10a)  Bj{t)uj{t)  —  bj{t)  £  ext  rge  Aj{t),  V<  G  JT.  j  =  1 . u, 

(4.106)  bj  =  -  Bj-iUj-i  -  .4j_iUj_,),  ;  =  1 . u. 

(ii)  A  differentiable  mapping  Xj  :  J  —*  .  j  =  0, . . .  .  i/  —  1.  is  a  solution  of  the  DAE 

(4.11)  Aj(t)Xj  4- Bj(t)xj  =  bj(t),  i£j,j=0 . -  1, 

if  and  only  if  x  j  has  the  form 

(4.12)  Xj=CjXj^t  +  Uj 
where  Xj+i  ;  Jf  — »  R'"-  is  a  differentiable  solution  of  the  DAE 

(4.13)  Aj+\{t)xj-^.i  +  Bj+l{t)xJ,^.l  =  bj,^.i{t),  t£j. 

(Hi)  A  differentiable  mapping  x  :  — ♦  R”  is  a  solution  of  (4.1)  if  and  only  if 

(4.14)  X  =  r^_ix„  +  iV-i, 
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where  r^_i  €  C'^(  JZ;  ,  K" ))  and  r„_i  €  C*  J7:  R")  are  given  by 

(4.15)  r,_,  =C7o  -  C  .1, 

(4.16)  Vi,-i  =  Co  ■  ■  •  +  Co  •  •  •  Cu-3^v-2  +  •  •  •  +  Cqui  +  Wo, 

and  Xi,  ■.  J  is  a  differentiable  solution  of  the  linear  ODE 

(4.17)  Ai,(t)x„  +  B^(t)x^  =  b^(t). 


Proof.  Parts  (i)  and  (ii)  follow  from  by  an  inductive  application  of  Theorem  4.1.  and  (iii) 
is  an  immediate  consequence  of  (ii).  □ 

The  essential  difference  between  (4.1)  and  (4.17)  is  that  the  latter  is  almost  an  explicit 
ODE  since  A^t)  €  GL{W''-'- )  except  perhaps  at  isolated  points.  Of  course,  the  case  when 
Auit)  €  GL{R.^‘'-'-)  for  every  f  €  JT”  is  of  special  importance.  For  this  case  we  have  the 
result: 

Corollary  4.2.  Under  the  same  hypotheses  and  with  the  notation  of  Corollary  4.1. 

(i)  the  condition 

(4.1S)  .4^0  €  Gi:(R’'‘'-‘), 

is  independent  of  the  reduction,  and 

(ii)  if  (4.18)  holds,  the  initial  value  problem 

(  A{t)x  4-  B{t)x  =  b(t)  in  J , 

(4.19)  \ 

[  z(^,)  =  z*,  €  v7,  z*  e  R", 

has  a  differentiable  solution  if  and  only  if 

(4.20)  z,  =  r^_i(<.)z^.  +  t;^_i(f.) 
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for  a  (unique)  x^.  €  .  Moreover,  if  (4.20)  holds,  the  differentiable  solution  x  of  (4.101 

is  unique  and  of  class 

Proof,  (i)  If  is  another  possible  choice  for  then  ~  (A^.B^)  by 

Theorem  3.1  (iii),  and  rank  A^{t)  =  rank  'it  ^  J  by  Theorem  3.1  (i). 

(ii)  If  (4.18)  holds,  then  (4.17)  is  an  explicit  linear  ODE.  and  hence  the  conclusion 
follows  from  Corollar\’  4.1  and  standard  existence  and  uniqueness  results  for  linear  ODEs 
-  which  in  particultir  implies  that  solutions  are  defined  in  all  of  J .  Note  that  uniqueness 
of  X ,,  in  (4.20)  follows  from  the  injectivity  of  r^_i(<, ).  □ 

Remark  4.1.  Condition  (4.20)  expresses  the  consistency  of  the  initial  value  r,  with  the 
DAE  at  □ 

For  b  =  0  in  (4.19)  we  may  choose  Uj  =  bj  =  0  in  Corollary  4.1  which  shows  that 
Vt/-i(l*)  =  0  in  (4.20).  Hence  it  follows  from  Corollary  4.2  (ii)  that  on  J7  the  homogeneous 
system 

(4.21)  Ax  +  Bx  =  0,  inj, 

has  ri,_i  linearly  independent  (and  analytic)  solutions  which  are  uniquely  determined  by 
their  value  at  any  point  €  J'.  We  now  show  that  the  hypotheses  of  Corollary  4.1  are 
necessary  for  the  existence  of  a  solution  for  cirbitrary  right-hand  side  in  (4.1),  and  that  the 
same  is  true  of  the  hypotheses  of  Corollary  4.2  (ii)  if  the  aforementioned  property  of  the 
solutions  of  (4.21)  is  to  be  preserved  as  well. 

Theorem  4.2.  (i)  In  order  for  the  DAE  (4.1)  with  analytic  coefficients  to  have  at  least 
one  differentiable  solution  in  J  for  every  b  €  C‘^(J7;R”).  it  is  necessar}'  that  (4.1)  be 
completely  reducible. 

(ii)  Let  (4.1)  be  completely  reducible  with  index  v  >  0.  If  the  homogeneous  problem 
(4.21)  has  r^_i  linearly  independent  differentiable  solutions  uniquely  determined  by  their 
value  at  any  given  point  t,  6  J.  then  At,{t)  G  GL{W''-^),  'it  ^  J  and  rank  Ay_i(t)  = 
'it  G  fj  if  V  ^  1. 

Proof,  (i)  It  suffices  to  prove  that  if  (4.1)  has  a  differentiable  solution  in  J  for  every 
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b  G  C‘*'(J7;R”),  then  the  pair  {A,B)  is  regular  in  and  the  reduced  DAE 

(4.22)  Ai(i)ii  +  B\{t)xi  =  bi{t),  ^  , 

has  at  least  one  differentiable  solution  in  for  every  6i  G  where,  of  course, 

r  =  majcjgj  rank  -4(^).  Indeed,  if  this  property  is  established,  it  may  evidently  be  used 
to  prove  inductively  the  complete  reducibility  of  (4.1). 

The  regularity  of  {A,B)  in  is  trivial  since  the  value  6(to)  <^3^^  arbitrarily  chosen 
for  any  given  aJid  existence  of  a  differentiable  solution  of  (4.1)  ensures  that  b{to)  G 
rge  .4(^o)  ©  B{to),  whence  rank  A{to)  ©  B{to)  =  n. 

Let  (i4i,5i)  be  a  reduction  of  {A,B)  so  that  (3.7)  holds  with  C  G  C‘^(j7’;  £(R’', R" )) 
and  D  G  C‘^(J7;  £(R”,R''))  satisfying  (3.5)  and  (3.6),  respectively.  For  any  6 1  G  C"{J';R’') 
the  existence  of  S  G  C‘^(j7’;  £(R'’,  R”))  such  that  Sit)  G  G£(R’’,  ext  rge  A{t)),  G  is 
guaranteed  by  Theorem  2.3  (i).  Clearly,  D(<)5(<)  G  GX(R’’),  G  J,  and  with  b{t)  = 
Sit)(D{t)Sit))~^bi{t),  'it  G  J,  we  obtain  b  G  C‘*'(j7’;R”)  such  that  b{t)  G  ext  rge  A{t)  and 
Dit)bit)  =  blit),  it  G  J. 

By  hypothesis,  the  DAE  (4.1)  has  a  differentiable  solution  x  for  this  6.  and,  in  addition, 
we  can  choose  uq  =  0  in  Theorem  4.1  (i).  Thus,  part(ii)  of  that  theorem  ensures  that  x 
has  the  form  x  =  Cxi  where  xi  is  a  differentiable  solution  of  (4.22)  in  J.  Thus.  (4.22) 
does  have  a  differentiable  solution,  as  claimed. 

(ii)  The  linearly  independent  solutions  of  (4.21)  can  be  viewed  as  the  columns  of 
an  n  X  differentiable  matrix  $(t)  satisfying  the  equation 

(4.23)  A(t)$  +  B(t)$  =  0,  t  G  J. 

Since  $(t)  is  uniquely  determined  (as  a  solution  of  (4.23))  by  its  value  $(t»)  at  any  t,  G  J . 
it  follows  that  $(t)  hcis  rank  r^_i  for  all  t  ^  J,  Since  we  may  choose  Uj  =  bj  =  0  when 
6  =  0  in  (4.1),  Corollary  4.1  provides  that  $(t)  =  r„_i(t)$„(t)  where  $i,  is  r„_]  x 
and  solves  the  system  .4^(t)4»^  +  B^(t)$^  =  0  in  Because  of  rank  $(t)  =  r^_i  we  also 
have  rank  $^(t)  =  and  therefore  $,/(t)  G  GI-(R''''"' ),  it  G  JT.  From  the  proof  of  part 
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(i),  the  existence  of  a  differentiable  solution  x  of  (4.1)  for  every  b  G  C''(J’;  IR")  implies 
the  existence  of  a  differentiable  solution  of  (4.17)  for  every  t  C~'(J:  R''"-'  ).  Clearlv 
Uu  —  is  differentiable  and  solves  the  system  A^{t)y^  =  f)^(t)  in  J  which,  in  turn, 

inplies  that  6^(f)  G  rge  A^{t).  Vf  G  J-  Since  b^,(t)  can  be  arbitrarily  chosen  for  given  t  ^  J 
this  shows  that  rank  .4^(0  =  r^_i  and  therefore  .4„(f)  G  Gi(R''‘'-“  ).  Vt  G  J . 

To  complete  the  proof,  we  have  to  show  that  rank  .4^_i(f)  =  r^_^.  Vt  G  J.  assuming,  of 
course,  u  >  1.  For  this,  it  is  obviously  not  restrictive  to  confine  attention  to  the  case  =  1 
and  hence  to  show  that  rank  A{t)  =  r{  =  tq),  V/  G  J-  The  necessary  conditions  already 
established  in  this  proof  enable  us  to  use  Corollary  4.2  with  6  =  0  which  shows  that  <f>  in 
(4.3)  is  not  only  differentiable  but  also  analytic  in  J.  By  Lemma  2.1  and  Theorem  2.1. 
with  .1/  =  we  obtain  L  G  C'^(  7:  0{n))  .such  that  U(t)  G  0(n)  is  a  linear  isomorphisms 

of  (rge  $(fo))'‘‘  onto  (rge  for  each  t  G  J.  Thus,  if  {ei . Cn-r)  denotes  a  basis  of 

(rge  the  matrix  ^(f ),  obtained  by  adding  the  column  vectors  L'(t)€,.  1  <  ?  <  n  — r. 

to  the  matrix  $(t).  has  full  rank  n  for  every'  i  €  J,  whence  4>  G  C‘^(^7;  GZ(R")). 

The  change  of  variables  y  =  ^~^x  transforms  (4.1)  into  the  DAE  A(/)4>(f  )i/-f  ( A(#)4>(/)4- 
Dit)^{t))y  =  6(f)  in  J,  which  therefore  has  at  least  one  differentiable  solution  for  every 
6  G  C“''(»7;E").  By  part  (i),  this  implies  that  the  pair  (.44..4l>  +  S4>)  is  regular.  On  the 
other  hand.  b\'  (4.23)  the  matrix  ,44>4-i?$  has  r  vanishing  columns.  Thus,  rank  ,4(f)4>( t)  + 
D{t)^{t)  <  n  —  r,  '7t  £  J .  But  then,  by  the  regularity,  w'e  must  have  rank  A{t)^{t)  >  r. 
Vf  G  J.  which  implies  that  rank  .4(f)  >  r  and  therefore  rank  .4(f)  =  r,  Vf  G  J .  □ 

We  note  here  that  some  of  our  arguments  in  the  proof  of  Part  (ii)  were  borrowed  from 
[CPS3]  and  [CS7].  Part  (ii)  of  Theorem  4.2  cannot  he  improved  to  guarantee  rank  .4^(f)  = 
7  =  0 . 1/  —  2.  for  all  f  G  ^7  as  the  following  2-dimensional  counterexample 

0  f  \ 

X  -P  X  =  6(  f ) 

0  0/ 

shows.  Tliis  D.A.E  has  index  2  and  the  unique  solution  (xi,X2)  =  (6)  —  tbi.  62).  In  this  case 
we  have  X]  =  0  and  (4.17)  is  the  trivial  equation  0  ==  0  in  R°  =  {0}. 

In  connection  with  Theorem  4.2.  the  setting  of  Corollary'  4.1  is  optimal  for  the  existence 
of  solutions  of  (4.1 )  for  arbitrary'  (  sufficiently  smooth)  right-hand  sides,  and  the  setting  of 
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Corollary  4.2  (ii)  is  optimal  if,  in  addition,  a  uniqueness  theory  similar  to  that  for  standard 
ODE  s  is  to  be  valid.  However,  Corollary  4.1  has  a  broader  range  of  applications  since  it 
also  allows  for  the  study  of  singularities  in  (4.1),  which  may  arise  if  the  matri.x  .4i,(C)  in 
(4.17)  becomes  singular  at  some  values  t  ^  J.  For  related  results  see.  for  instance.  [H64] 
or  [F65]. 

At  this  stage,  it  is  important  to  note  that  the  reduction  procedure  of  Section  3  for 
analytic  pairs  (A.B)  does  not  require  regularity  of  (A,B)  and  indeed  can  also  be  performed 
under  the  weaker  assumption  that  maxt^j  rank  .4(f)  0  B(f)  =  n.  In  particular,  the 
relation  (3.4)  remains  true  cifter  enlarging  the  set  S  so  as  to  incorporate  the  (isolated) 
points  where  rank  .4(f)  ©  B{t)  <  n.  .A,U  the  results  of  Section  3  remain  valid  with  the 
Scune  proofs  if  complete  regularity  is  now’  replaced  by  the  more  general  assumption  that 
maxjgj  rank  Aj{t)  ©  Bj{t)  =  rj_i.  (r_i  =  n).  But  under  these  weaker  assumptions,  the 
D.4E  (4.1)  need  no  longer  be  reducible,  for  the  conclusion  of  Theorem  4.1  (i)  does  not 
remain  true  for  arbitrary  right-hand  sides  b.  As  a  consequence,  for  the  \'alidity  of  Theorem 
4.1  (ii)  it  can  no  longer  be  proved,  but  must  be  assumed,  that  there  exists  a  uq  €  C^{J\ R") 
for  w’hich  B(t)uo(t)  —  b{t)  €  ext  rge  .4(f).  .4  similar  assumption  must  be  made  about  the 
bj's  in  Corollary  4.1  to  arrive  at  the  reduced  problem  (4.17).  In  other  words,  the  reduction 
of  (4.1)  remains  possible  only  for  restricted  choices  of  right-hand  sides,  independent  of  *^he 
smoothness.  To  give  an  idea  of  the  nature  of  the  required  restrictions,  w’e  mention  w’ithout 
proof  that  if  6  G  C^'lJT^iR"),  then  Theorem  4.1  (i)  remains  v-alid  if  the  zeroes  of  the  vector 
Bit)^ Adj[A{t)Ait)^  —  B{t)B{t)^]b{t)  -  where  ''Adj''  stands  for  the  “.4djugate";  that  is. the 
transpose  of  the  “matrix  of  cofactors"  -  ail  have  order  greater  than  or  equal  to  the  order 
of  the  zeroes  of  det(.4(f).4(f )^  +  B{t)B{t)^). 

Finally,  w’e  note  that  a  generalization  of  the  reduction  procedure  of  Section  3  for  analytic 
pairs  [A.B)  is  always  possible  w’ithout  any  additional  assumptions.  However,  then  the 
reduced  pairs  {Aj.Bj)  axe  such  that  Aj,Bj  €  C“'(77;  .R'"^"’ ))  with  rij  >  in 

other  words.  Aj  and  Bj  need  no  longer  be  square  matrices.  Once  again,  the  reduction 
stabilizes  after  a  finite  number  n  >  0  of  steps,  up  to  a  straightforw’ard  generalization  of 
the  equivalence  of  Section  3.  But.  in  the  case  when  >  Cj  . j  for  at  least  one  index  j. 
it  can  be  shown  that  local  uniqueness  for  the  solutions  of  the  initial  value  problem  (4.19) 
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is  not  true.  This  is  due  to  the  fact  that,  when  possible,  the  reduction  to  the  form  (4.17) 
leads  to  a  system  with  more  unknowns  than  equations  and  a  surjective  .4^(t)  except  at 
isolated  points.  On  the  other  hand,  the  condition  Vj  >  0.  amounts  exactly  to 

maxfgjrcink  Aj{t)  0  Bj{t)  =  Tj-i,  Vj  >  0.  the  case  discussed  above.  For  some  further 
comments  along  this  line  see  also  Section  7. 

5.  Global  reduction  of  differential-algebraic  equations:  generalized  solutions. 

Let  T>{J)  denote  the  space  of  infinitely  differentiable  real  valued  functions  with  compact 
support  in  J  and  V'{J)  its  dual,  the  space  of  distributions  in  J .  For  x  =  ( j.)  ^  [V{J)Y 
and  ^  e  {V[J)Y,  we  set 


n 

=  i: 

t=i 

where,  on  the  right.  (•,  •)  denotes  the  duzdity  pairing  between  ViJ)  and  V{J).  For 
M  €  C~(J;£(R”.Il"*))  and  x  €  the  product  Mx  e  (P'(J))"'  is  defined  by 

We  consider  now  the  generafized  DAE 
(5.1)  A{t)i  +  B{t)x  =  b,  6G(r'(J))". 

where,  as  before,  the  coefficients  A  and  B  are  analytic  in  J.  A  solution  of  (5.1)  is  any 
X  €  ('D'(J7’))"  such  that  (.4x  +  Bx,^p)  =  (6, y?).  for  all  G  [T>{J))'^.  where,  of  course,  x 
denotes  the  derivative  of  x  in  the  sense  of  distributions. 

For  the  extension  of  the  reduction  of  Section  4  to  (5.1)  we  need  a  preliminaiy  lemma. 

Lemma  5.1.  Suppose  that  the  pair  (A.B)  is  reguJai-  and  set  r  =  maxtejrank  A{t).  Let 
P  G  £(K"))  be  such  that  P{t)  is  a  projection  onto  ext  rge  .4(f).  Vf  G  J,  and  denote 

by  n(f)  the  orthogonal  projection  onto  ker  Q(t)B(t}.  Vf  G  57’,  where  Q  =  I  —  P .  Then: 

(i)  For  every  --p  G  there  exists  a  0  G  such  that  [I  —  =  B^Q^v. 
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fh)  Let  C  e  £(R'',R”)}  satisfy  (3.5).  Then,  for  x  €  {V{J))’',  there  exists  a 

xi  G  {V{J)Y  such  that  x  =  Cxj  if  and  only  if  (I  —  n)i  =  0. 

(Hi)  Let  D  G  C‘*’(  J”;  £(R”,1R’'))  satisfy  (3.6).  Then,  for  w  G  (T>iJ)Y.  there  exists  a 
G  €  {V(J)Y  such  that 

Proof,  (i)  In  line  with  Remark  3.1  we  have  dimker  Q(t)B(t)  =  r  for  t  G  JT"  and  hence  it 
follows  that  rge  B{t)^Q(t)'^  =  [ker  Q( 0-5(0] =  rge  (7  —  11(0)  's  (n-r)-dimensional.  and. 
because  of  rank  Q(t)^  =  rank  Q(t)  n  ~r.  that  B(t)^  G  GL{rge  Q{t)^ ,  [ker  Q(05(0]"*’)- 
By  Theorem  2.2  (i)  and  (hi),  there  axe  mappings  S  G  C‘^(J7;  £(R’’'''’,  R" ))  and  V'  G 
C^(  J;£(R",R'*-0)  such  that 

5(i)  S  GI(R"-',rge  Q{tf  ),  V-(<)|,„  €  GI(rge  BitfQitf.W-') 

and  kerF(t)  =  kerg(05(0,  Vt  G  J.  This  implies  that  VB'^S  G  C'^{J\GL(W^-^)). 
Thus,  for  G  (T*(*7))")  the  choice  if  =  S(V’B^S)~^  provides  that  and 

VB^Y  =  V<f\  that  is,  (I  —  Tl)(B^Q^Y  —  v’)  =  0.  But  since  (7  —  ]T)B^Q^  =  B^Q^.  this 
means  that  B^Q^Y  =  —  n)v?. 

(ii)  Suppose  first  that  x  =  Cxi-  Since,  again  in  view  of  Remark  3.1,  kerC(0^  = 
[rge  C(0]'*'  =  [ker  (5(05(01  we  have  C^(7  —  II)  =  0  and  therefore.  {x,(7  —  IT);,?)  = 
(x, .  C^(7  -  n)(,?)  =  0.  G  (I>(  J))";  that  is,  (7  -  n)x  =  0. 

Conversely,  suppose  that  (7  —  n)x  =  0  and  set  Xj  =  (C^C)~^  C^x.  Then,  for  G 
iV{J)Y.  we  have  {Cx,,^)  =  {x,C{C^C)-^C’^^).  But,  C(0(C(0^C(0)"'C'(0^  is  the 
orthogonal  projection  n(t)  onto  rge  C{t)  =  ker<5(t)5(t)  (see  once  again  Remark  3.1). 
Thus,  because  (x,(7  —  n)^?)  =  0  by  hypothesis,  we  find  that  (Cxi,(,?)  =  (x,!!:,;)  =  (x.^) 
which  gives  x  =  Cx  i . 

(iii)  If  we  can  show  that  P(t)^  G  GL(rge  D{t)^,  rge  P(O^),  V<  G  J.  the  method  used  in 

the  proof  of  (i)  above  will  provide  the  desired  result.  Since  rank  D(t)^  =  rank  D{t)  =  r. 
and  rank  5(0^  =  rank  P(t)  =  r.Yt  e  J.  it  suffices  to  show  that  P{t)'^ u  =  0  for  some 
u  G  rge  D(t)^  implies  that  u  =  0.  Such  a  u  is  an  element  of  [ker  7)(f)]-^ n  [ext  rge  -df/)]-*-  = 
[ker  D(t)  +  ext  rge  .4(0]‘^-  But,  since  kerD(0  C  ext  rge  ,4(0  =  {0},  we  have  ker  D(t)  9 
ext  rge  A(t)  =  R".  so  that,  indeed,  u  must  be  0.  □ 
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We  are  now  in  a  position  to  prove  the  “distribution"  version  of  Theorem  4.1. 


Theorem  5.1.  Suppose  that  the  DAE  (5.1)  nith  analytic  coefEdents  and  r  =  max(gj.4(t) 
is  reducible,  and  let  (,4i  ,Bi)  be  any  reduction  (3.7)  of  the  pair  {A.  B)  in  J  with  mappings 
C  G  £(R'',R")),  and  D  G  £(R",  R'"))  satisfying  (3.5)  and  (3.6),  respectively 

Then: 

(i)  There  exists  a  uq  G  {V{J)Y  such  that 

(5.2)  Q(Buo-b)  =  0, 

where  Q  =  I  —  P  and  P  G  £(R"))  is  such  that  P(t)  is  a  projection  onto  ext  rge  .4(/). 

'it  ^  J .  In  particular,  (5.2)  holds  for 

(5.3)  uo  =  B'^iAA'^  +  BB'^r^b. 

(ii)  A  distribution  x  G  (D'{J)Y  solves  (5.1)  if  and  only  if  for  any  uq  satisfying  (5.2) 

(5.4)  X  =  Cxi  +  Uo, 

where  xj  G  solves  the  DAE 

(5.5)  i4iXi  +  BjXi  =  bi, 
where  b\  is  given  by 

(5.6)  bi  =  D(b  —  Buo  —  Aiio). 

Proof,  (i)  Let  uq  be  defined  by  (5.4)  and  set  uq  =  A^{AA^  +  BB^)~^b.  Then,  .4t’o  +  Buo  = 
b  and  for  G  (T>(J)Y  we  have  {Q{Buo  -  b),ip)  =  {Buq  -  b,Q'^ip)  =  -{Avo.Q'^g>)  = 
—  {I'o.A^Q^'g))  =  0  since  A^Q^  =  [QA]^  =  0. 

(ii)  Let  X  G  {D{J)Y  be  a  solution  of  (5.1),  so  that  for  any  ip  G  (D{J)Y  we  have 
{Q{Bx  —  b),ip)  =  {QiAx  -f  Bx  —  b),ip)  =  {Ax  +  Bx  —  b,Qip)  =  {O.Qip)  =  0  and  therefore 
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Q{Bx  —  6)  =  0  which,  because  of  Q{Buq  —  6)  =  0.  implies  that  QB(x  —  uq)  =  0.  Thus 
(x  —  Uq,  B^Q^K')  =  0.  for  edl  i/j  €  and,  hence,  by  Lemma  5.1  (i),  (x  —  uq,{I  — 

n)>?)  =  0,  Vi,?  €  (P(v7))";  that  is,  (/  —  n)(x  —  uq)  =  0.  Thus,  Lemma  5.1  (ii)  ensures 
that  X  =  Cxi  +  Uq  for  some  Xi  G  {'D'{J)Y.  Substituting  this  into  (5.1),  we  find  that 
TCxi  +  {BC  +  /lC)xi  —  h  —  Buq  —  i4iio,  and  multiplying  both  sides  by  D  we  obtain  that 
Xi  solves  (5.5). 

Conversely,  let  xi  G  {T>'{J)Y  solve  the  DAE  (5.5)  and  set  x  =  Cxi  +  uq,  whence 
D{Ax  +  Bx  —  b)  =  0.  By  construction  of  x,  we  have  Q{Bx  —  b)  —  Q{BCxi  -)-  Buq  —  6)  =  0 
since  QBC  =  0  (see  Remark  3.1)  and  Q{Buq  —  fc)  =  0.  Since  QA  =  0.  this  implies  that 

(5.7)  Q{Ax Bx  —  b)  =  Q, 

and  therefore  DP{Ax  +  Bx  —  b)  =  D(Ax  +  Bx  —  b)  =  0  which  is  equivalent  with 

(5.8)  (Ai +  Bx -b,  =0,  V<pe(V(J}Y- 

Let  G  (’D(J7’))",  so  that  tp  —  P'^rp  +  Q'^tp.  By  Lemma  5.1  (iii),  there  exists  a  G  ['P{J)Y 
such  that  P^xp  =  P^D^tp.  Thus,  xp  =  P^D^xp  +  Q^xp,  and  (Ax  +  Bx  —  b,xi')  =  {Ax  + 
Bx  —  b,  P'^ D^xp)  +  (Ax  4-  Bx  —  b,Q'^xp)  =  0  by  (5.7)  and  (5.8).  This  shows  that  x  solves 
(5.1).  □ 

By  recursive  application  of  Theorem  5.1,  we  see  that  when  the  DAE  (5.1)  is  completely 
reducible,  its  solutions  are  of  the  form  x  =  r„_iXi,  +  where  x^  G  {'P'{J)Y‘'~'  solves 

(5.9)  A„x„  +  B„x„  =  bj,, 

and  A^{t)  is  invertible  for  all  but  isolated  values  of  t  G  57.  Here,  A^,  B^,  r„_i  and  are 
obtained  as  in  Corollary  4.1  with  the  obvious  modifications  that  now  Uj,bj  G  i'D'iP7)Y^~' 
(hence  tv-i  ^  {P'{J)Y)  and  that,  instead  of  (4.10a),  Uj  is  required  to  satisfy  the  condition 
Qj{BjUj  —  bj)  =  0,  where  Qj  =  I  —  Pj  and  Pj  G  C{W’-'^ ))  is  such  that  Pj{t)  is  a 

projection  onto  ext  rge  Aj(^),  Vt  G  57. 
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If  A^(t)  is  invertible  for  every  t  ^  J ,  then  (5.9)  is  simply  the  explicit  ODE  with  analytic 
coefficients 

whose  solutions  form  an  rj,_i -dimensional  affine  subspace  of  (see  e.g.  Schwartz 

[S66]).  The  solutions  of  (5.1)  thus  form  an  r^_i -dimensional  affine  subspace  of  (X>'(J7))". 
Of  course,  in  this  framework,  initial  value  problems  make  no  sense,  in  general,  since  dis¬ 
tributions  need  not  have  pointwise  values.  However,  Corollary  4.2  can  be  extended  for 
special  choices  of  h.  as,  for  example,  b  G  foi'  some  k  >  v.  In  this  case  we  note, 

without  proof,  that  Corollary  4.2  has  an  obvious  generalization  due  to  the  fact  that  the 
solutions  of  (4.1)  eire  in  {J)  C  C^{J)  and  hence  have  numerical  values  at  all  the 

points  of  J . 

Initial  value  problems  also  malce  sense  when  the  right-hand  side  of  (5.1)  is  a  discontin¬ 
uous  function,  provided  that  the  initial  value  is  not  chosen  at  a  point  of  discontinuity.  It 
turns  out  that  this  result  yields  rather  complete  answers  to  the  problem  of  solving  (4.1) 
with  “inconsistent”  initizd  condition;  that  is,  with  an  initial  condition  i,  that  does  not 
satisfy  the  requirement  (4.20).  This  problem  has  been  discussed  at  some  length  in  the 
constant  coefficient  case  because  of  its  relevance  in  problems  from  various  applications. 
Details  will  be  presented  elsewhere. 

Remark  5.1  At  the  end  of  Section  4,  we  observed  that  the  reduction  procedure  for 
the  DAE  (4.1)  could  be  extended  to  any  system  with  analytic  coefficients,  provided  that 
suitable  limitations  are  placed  on  the  right-hand  sides  bj.  It  is  interesting  to  note  that 
this  is  no  longer  true  if  the  problem  is  understood  in  the  sense  of  distributions.  Indeed, 
if  the  pair  {A,B)  is  not  regular,  part  (i)  of  Lemma  5.1  fails  to  hold  even  in  the  czise 
when  maxegjrank  A(t)  @  B{t)  =  v  ond  hence  when  kev  Q{t)B{t)  must  be  replaced  by 
rest  ker  Q{t)Bit).  Indeed,  at  a  point  to  €  ,7  where  rest  keTQ{to)B{to)  C  ker(5(to)-B(to). 
we  have  rge  B(to)^Q(to)^  £  rge  (/  —  n)(to),  and  if  9  €  {V{J)Y  is  chosen  such  that 
(/  —  n)(to)9(to)  7  B{to)^Q{to)^ ,  it  is  clear  that  no  ^  G  {'D{J)Y  exists  such  that 

(/  —  11)9  =  B^Q^v.  But,  without  Lemma  5.1  (i),  it  is  no  longer  possible  to  ascertain 
that  every  solution  of  (5.1)  must  have  the  form  (5.5),  even  if  existence  of  uq,  as  in  (5.2), 
is  assumed.  The  converse,  however,  remains  true;  that  is,  if  xi  solves  (5.6),  then  x  given 
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by  (5.5)  solves  (5.10).  This  discrepancy  may  be  explained  by  the  fact  that  the  crucial 
implication  (ii)  =>  (i)  in  Lemma  3.2  relies  heavily  upon  a  continuity  argument  for  K 
(meaningless  when  K  =  x  —  uq  is  a  distribution)  if  rest  ker  Q{t)B{t)  and  kcTQ{t)D(t)  do 
not  coincide  for  all  t  ^  J .  □ 

6.  Boundary  value  problems. 

The  global  reduction  result  of  Corollary  4.1  (iii)  is  especially  convenient  for  handling 
boundary  value  problems.  Indeed,  suppose  that  we  associate  the  completely  reducible 
DAE  (4.1)  with  index  v  with  a  condition  of  the  form 

(6.1)  5(a,6,r(a),r(6),i(a),i(6))  =  0, 

where  g  :  x  JT"  x  R"  x  R"  x  R"  x  R"  — ♦  R"*,  and  a, 6  6  J-  By  Corollary  4.1  (hi)  (cmd 
using  the  notation  of  that  corollary)  we  know  that  r  :  — »  R"  solves  (4.1)  if  and  only  if 

(6.2)  X  =  r^^iXu  +  Vu-u 
where  x^  :  J  R''*--'  solves  the  ODE 

(6.3)  x^  +  B^{t)xi,  = 

assuming,  of  course.  .4^(0  €  GL{W‘’-^),  Vt  6  »7-  Thus,  by  substituting  (6.2)  into  (6.1). 
we  see  that  x  :  >  R"  solves  (4.1)  and  satisfies  the  boundary  condition  (6.1)  if  and  only 

if  X  has  the  form  (6.2)  with  a  solution  x„  of  the  ODE  (6.3)  that  satisfies  the  boundary 
condition 

(6.4)  g^(a,6,x^(a),x^(6),x„(a),x^(6))  =  0, 
where  :  J  x  J  x  R’’--*  x  R'""-*  x  R'’*'-'  x  R’’*'-'  —*■  R"*  is  given  by 

<7^(#,r,x^,p^,,^^.r?^)  =  <7(t.r,r„_i(0x^  +  t;^_,(t),  r^_i(0Pi/  +  f^_,(f)x^  +  tV-i(0. 

(6.5)  r„_i(r)^^  +  u„-i(r),r„_i(r)r7^  +  f„_i(r)^^  +  tv-i(r)). 
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But  in  fact,  since  by  (6.3),  x„(a)  eind  i„(6)  are  given  in  terms  of  a,  x^(a)  and  6,  x^(b), 
respectively,  the  boundaiy"  condition  (6.4)  takes  the  form 

(6.6)  h^{a,b,x^(a),x^(b))  =  0, 

where  :  J  x  J  x  R'"*'-'  x  R''*'-^  —*  R"*  is  obtained  by  substituting 

Pk  =  Au{t)~'^[b^it)  -  T]u  =  A^(r)"^[6^(r)  - 

into  the  arguments  of  g„  in  (6.5). 

The  problem  (6.3)  with  the  boundzu’y  condition  (6.6)  is  a  boimdairy  value  problem  for  an 
explicit  linear  ODE  which  therefore  can  be  handled  by  the  standard  theory.  Furthermore, 
if  g  =  g{t,T,x,p,^^Tj)  in  (6.1)  is  an  affine  function  of  (x,p,^,g),  then  h^,  = 
in  (6.6)  is  an  affine  function  of  (x„,(u)- 

The  exact  same  considerations  apply  to  the  distribution  case  if  6  6  ^^lo’c  (^7)-  Indeed,  if 
so,  we  may  guarantee  that  the  sequence  {bj)  of  Corollary  4.1  satisfies  bj  €  for 

0  <  i  <1^,  whence  r„  €  (J)  C  C^{J)  has  a  numerical  value  at  each  point  of  J . 

7.  The  non-analytic  case. 

For  .4.  B  G  C°°(j7;  T(R")),  the  regularity  condition  (3.2)  is  no  longer  sufficient  to  repeat 
the  arguments  of  Section  3  and  define  reductions  (.4i,Bi)  of  {A,B).  This  is  due  to  the 
fact  that  the  concepts  of  extended  range  and  restricted  null-space  do  not  necessarily  make 
sense,  and  that  rank  .4(t)  need  no  longer  equal  its  ma^cimum  value  in  J  at  all  but  isolated 
points.  However,  if  we  strengthen  the  regularity  condition  as  follows 

(7.1)  rank  ^(t)  0  B{t)  =  n,  rank  A{t)  =  r,  V<  G  .7, 

then  reductions  (Ai,Bi)  of  (A.B)  can  be  obtained  by  the  method  of  Section  3.  To  see 
this,  it  suffices  to  replace  “ext  rge”  and  “rest  ker”  by  “rge"’  and  “ker’\  respectively,  and 
to  use  Theorem  2.2  in  lieu  of  Theorem  2.3.  If  (7.1)  again  holds  for  the  pair  (Ai,Bi) 
with  n  replaced  by  r,  and  r  replaced  by  a  new  integer  rj  <  r,  then  we  obtciin  a  new 
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pair  (.42,i?2)  by  reduction  of  (.4i,i7i),  and  so  on.  The  proof  of  Theorem  3.1  can  be 
repeated  verbatim,  provided  that  is  replaced  by  everywhere,  including  in  the 

equivalence  definition  of  pairs  (Definition  3.1).  Likewise,  Theorem  3.2  has  an  obvious 
extension  if  complete  regularity  of  the  pair  (A.B)  is  replaced  by  the  stronger  requirement 
that  for  j  >  0  the  pair  {Aj,Bj)  satisfies  the  appropriate  analog  of  condition  (7.1).  In 
particular,  an  index  >  0  can  be  assigned  to  the  pair  (A.B)  in  this  case. 

In  line  with  this,  it  should  be  evident  how  the  results  of  Section  4  up  to  and  including 
Corollary  4.2  and  those  of  Sections  5  and  6  can  be  extended  when,  with  .4o  —  A.  Bq  =  B . 
we  have 

(7.2)  Tai]k  Aj(t)  Q  Bj{t)  =  Tj-i,  Vf  €  i/,  Vj  >  0, 

(7.3)  rank  Aj(t)  =  rj,  \/t  €  »7,  Vj  >  0, 

where  r_i  =  n  >  Tq  >  ■  ■  ■  >  Tj  >  and  for  each  j  >  1,  {Aj^Bj)  denotes  a  reduction 

of  {Aj-i,Bj-i).  On  the  other  hand,  Theorem  4.2  has  no  interesting  generalization  to  this 
case,  and  most  of  the  discussion  concluding  Section  4  becomes  immaterial  with  one  notable 
exception;  The  reduction  of  peurs  {A,B)  remains  possible  with  Aj,  Bj  being  rij  x  rj_i. 
if  (”-3)  is  unchanged,  and  (7.2)  is  replaced  by  the  weaker  condition 

(7.4)  rank  Aj[t)  ©  Bj{t)  =  pj,  Vf  €  J.  r_,  <  Pj  <n,  Vj  >  0. 

Once  again  there  is  a  smallest  integer  i/  >  0  such  that  and  -4^(0  is  onto  for 

all  f  G  As  in  Section  4,  if  Hj  >  rj_i  for  at  least  one  j,  local  uniqueness  of  the  solutions 
of  (4.19)  is  lost.  Furthermore,  because  of  (7.4)  it  is  etisy  to  characterize  those  right-hand 
sides  h  for  which  the  sequences  Uj,  bj  are  defined,  and  hence  for  which  the  reduction  of 
(4.1)  to  the  form  (4.17)  is  possible.  In  fact,  it  is  both  necessary  and  sufficient  that  (with 
ho  =  h) 

(7.5)  hjit)  e  rge  .4,(O0B,(O  =  rge  {Aj{t)Aj{tf  +  B,it)Bj(tf),  t  e  J .  0  <  j  <  -  1. 

where  the  surjectivity  of  Ai,{t)  ensures  that  this  condition  holds  for  j  =  v.  If  so.  there  exists 
a  Wj  as  smooth  as  bj  such  that  bj  =  {Aj{t)Aj(t)^  +  Bj{t)Bj{t)^)wj,  and  we  may  choose 
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U]  =  BJwj.  Obviously,  the  condition  (7.5)  is  also  necessary  for  existence  of  a  solution.  If 
6  is  a  distribution,  the  condition  (7.5)  is  replaced  by  bj  =  (Aj(t)Aj(t)^  +  )u'j 

where  Wj  6  (P'(R))";  we  shall  not  enter  into  the  details. 

Similar  conclusions  regarding  this  “Fredholm  Alternative"  for  undetermined  system 
(tIj  >  rj_i  for  at  least  one  j)  axe  obtained  in  [KuMe92]  assuming  three,  instead  of  our 
two,  constant  ramk  conditions  at  each  step  of  their  reduction. 

Remark  7.1.  When  .4  and  B  axe  analytic  and  the  DAE  (4.1)  has  index  1,  it  follows  from 
Theorem  4.2  that  both  conditions  in  (7.1)  as  well  as  rank  .4](t)  =  r,  Vt  G  are  necessary 
for  the  validity  of  an  existence  and  uniqueness  theorj'  similar  to  that  for  ODE's.  Thus, 
the  hypotheses  of  this  section  are  identical  to  those  made  in  the  analytic  case  for  index  1 
problems  “without  singularities".  This  is  no  longer  true  for  higher  index  problems.  □ 

Extending  the  results  about  classical  solutions  from  A,R  G  C°°(  £(R" ))  to  A,  B  G 

C'"(  £(R")),  m  >  n,  requires  nothing  more  than  replacing  by  “C"”’  everj’where. 

and  u  <  k  <  oo  by  u  <  k  <  m  in  Theorem  4.1  as  well  as  Corollaries  4.1  and  4.2.  In 
fact,  even  C"  regularity  is  mostly  superfluous  tmd  only  C''  regulzirity  is  needed  where  u 
is  the  index  of  {A.B).  But  since  v  (tind  even  the  existence  of  i/)  is  not  known  a  priori, 
the  hypothesis  A,B  G  C‘'(.7;  £(R"))  is  ambiguous  and  should  be  replaced  by  “.4.  B  are 
sufficiently  smooth". 

A  little  more  care  must  be  exercised  regarding  generalized  solutions.  Recall  that,  given 
an  integer  k  >  0,  n  distribution  i  G  {V'{J)Y  is  said  to  be  of  order  k  if  (x.  •)  extends 
by  continuity  to  elements  of  (X>*'(J7))",  where  is  the  space  of  real-valued.  A;-times 

continuously  differentiable  functions  in  J  with  compact  support  (see  [S66]  for  details). 
Distributions  of  order  k  can  be  multiplied  by  (matrix)  functions  ii  C  >  k  but  not  if 
[  <  k.  .4s  a  result,  if  .4.5  G  C”’(J7:  £(R"))  htis  index  u  in  the  sense  of  this  section  (so  that 
necessarily  m>  v).  then  Au  and  B„  are  of  cla.ss  C'"”",  and  is  of  order  k  u.  In  order 
to  transform  (5.10)  into  an  explicit  ODE  by  inverting  -4„,  we  must  be  able  to  multiply  b^ 
by  .4(7*  -  Since  .4j7'  is  C^~‘' ,  this  requires  m  —  u  >  k  A  v  and  hence  k  <  m  —  2v.  In  other 
words,  in  general,  the  reduction  of  Section  5  can  be  performed  only  if  m  >  2v,  and  only 
when  the  right-hand  side  6  is  a  distribution  of  order  m  —  2v  or  less.  On  the  other  hand,  if 
b  G  IT|^7^(77)  it  is  easily  seen  that  m  >  suffices  as  in  the  case  of  classical  solutions. 
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Problems  with  index  1  have  been  studied  extensively,  at  least  w'hen  no  singularities 
are  present  and  when  attention  is  confined  to  classical  solutions.  Since  the  hypotheses  for 
index  1  used  in  this  section  prohibit  existence  of  singularities,  it  may  be  of  interest  to  check 
how  they  relate  to  others  made  in  the  literature.  Throughout  this  discussion  the  stronger 
hypotheses  of  this  section  are  used  under  which  '"ext  rge"  and  “rest  ker "  become  "rge  " 
and  “ker  ”,  respectively. 

Theorem  7.1.  For  A.B  €  C'(»7;£(R"))  the  pair  (A,B)  has  index  1  in  J  in  the  sense  of 
this  section  if  and  only  if 
(i)  rank  Ait  j  =  r,  'it  ^  J 

and  one  of  the  following  four  equivalent  conditions  hold: 

(ii-a)  {i  G  ker .4(f),  B(t)x  £  rge  >l(f)}  =>  a:  =  0,  it  £  J , 

(ii-b)  A  +  QB  £  C\J\  GL(R^)),  where  Q  =  I  -  P  and  P  £  C^(J7;£(R”))  is  such  that 
P{t)  is  a  projection  onto  rge  .4(f),  it  £  J . 

(ii-c)  A  +  QiB  +  .4)  €  CHJ;  GI(R”))  with  Q  as  in  (il-b). 

(ii-d)  A  +  BK  £  C'(J7’;  GT(R”)),  where  K  €  C^( JT;  £(R’’))  is  such  that  K(t)  is  a 
projection  onto  ker  .4(f),  it  £  J . 

Proof,  (a)  Suppose  that  {A,B)  has  index  1  in  ff,  so  that  (i)  holds.  Using  the  previous 
notation,  we  have  .4i(f)  £  GL{W),  it  £  J whence  by  (3.6)  and  (3.7),  A(t)Cit)  has 
rank  r  for  t  £  J .  By  (i)  this  occurs  if  and  only  if  rge  C(f)  D  ker -4(f)  =  {0}.  But 
rge  C(f)  =  keTQ{t)B(t)  by  (3.5),  and  hence  we  have  ker  Q(f)F(f)nker .4(f)  =  {0}.  it  £  J . 
Since  Q(t)  projects  onto  a  complement  of  rge  i4(f),  it  follows  that  x  £  ker  Q(f)B(f  inker  .4(f) 
if  and  only  if  x  G  ker  .4(f)  and  B(t)x  G  rge  Ait)  which  proves  (ii-a).  Conversely,  assume 
that  (i)  and  (ii-a)  hold.  Condition  (ii-a)  and  the  second  part  of  (7.1)  coincide.  By  (i)  and  (ii- 
a),foreachf  G  if  the  null-space  of  the  mapping  (p,i)  G  R”xR"  i— >  (j4(f)x.  .4(f)p-t-i?(f  )x)  G 
R"  X  R"  has  dimension  dim  ker  .4(f)  =  n  —  r.  Hence,  its  rank  equals  p  =  n  +  r.  But 
clearly  p  <  rank  .4(f)  -(-  rank  .4(f)  ©  Bit)  =  r  -t-  rank  Ait)  ©  Bit).  This  implies  that 
rank  .4(f)  ©  F(f)  =  n.  it  £  J  which  proves  the  first  part  of  (7.1)  and  allows  us  to  speak 
now  of  a  reduction  (.4i,5i)  of  iA,B)  in  J.  Then,  reversing  the  steps  in  the  first  part  of 
the  proof,  we  arrive  at  rank  .4i(f)  =  r,  Vf  G  if  ;  that  is,  the  index  of  (.4, B)  in  J  is  1. 


34 


(b)  By  Lemma  2.1.  condition  (i)  ensures  the  existence  of  P  €  C^(v7:£(R"))  such  that 
P(t)  is  a  projection  onto  rge  -4(0,  V#  6  J7,  (actually  it  is  the  orthogonal  projection).  With 
Q  =  I  —  P  it  IS  obvious  that  (ii-a)  is  equivalent  with  (ii-b). 

(c)  In  order  to  show  that  conditions  (ii-b)  and  (ii-c)  are  equivalent,  suppose  first  that 

(ii-b)  holds  and  let  u  €  R"  be  such  that  A{t)u  -f  A{t))u  =  0  wliich.  obviously, 

requires  that  u  G  ker-4(0-  Now,  using  QA  =  0  we  see  that  QA  =  —QA.  and  hence  that 
Q{t)B{t)u  =  0.  Thus.  A{t)u  -t-  Q{t)B{t)u  =  0,  whence  u  =  0  and  (ii-c)  holds.  Conversely, 
if  (ii-c)  is  valid  then  for  u  G  ker(j4(0  +  Qii)B{t))  we  have  u  G  ker-4(0  H  ker  Q{t)B{  t) 
and,  because  of  u  G  ker-4(0,  it  follows  that  u  G  keTQ{t)A{t)  =  kerQ(0^(t)-  Thus. 
u  G  ker  (-4(0  +  Qi^)B{t)  -|-  ^(0-4(0)  and  (ii-c)  implies  that  u  =  0  which  proves  the 
validity  of  (ii-b). 

(d)  For  the  equivalence  of  (ii-b)  and  (ii-d)  assume  first  that  (ii-b)  is  valid  and  let  u  G  R" 

be  such  that  A{t)u  -|-  B{t)K{t)u  =  0  and  hence  Q{t)B{t)K{t)u  =  0.  Then  we  have 
[-4(0  +  Q{t)B{t)]K{f)u  =  0,  since  -4A'  =  0,  and  therefore  K{t)u  =  0  and  A(t)u  —  0: 
that  is,  u  G  ker-4(0-  But  then  u  =  K{t)u  =  0  shows  that  (ii-d)  holds.  Conversely, 
suppose  that  (ii-d)  holds  so  that  A{t)'^  -I-  K{t)^B{t)^  G  GI.(R").  Since  K{t)^  projects 
onto  a  complement  of  rge  .4(0^,  it  plays  the  szune  role  as  Q(t)  above,  but  with  -4(0> 
Bit)  replaced  by  -4(f)^,  P(0^,  respectively.  Likewise,  Q{t)^  plays  the  role  of  K{t)  in 
this  case.  Thus,  from  the  first  part,  we  have  A{t)^  +  B{t)^Q{t)^  G  GL(R"),  and  hence 
-4(0  +  G  GL(R")  by  transposition.  This  shows  that  (ii-b)  holds.  □ 

Conditions  (i)  and  (ii-a)  are  exactly  the  index  1  conditions  in  the  sense  of  [RRh92] 
for  the  D-A.E  (1.1).  The  simplest  way  to  see  this  is  to  use  the  remark  in  [RRh92]  that 
“index  1”  in  the  sense  of  that  paper  is  the  same  as  “index  T'  in  the  sense  of  [RRh91],  and 
that  the  conditions  given  in  the  latter  paper  are  equally  valid  without  the  second-order 
derivative  terms.  With  this  simplification,  the  constant  rank  condition  and  condition  (2.22) 
of  [RRh91]  -  which  are  necessary  and  sufficient  for  index  1  -  become  rank  DpF(t.  x,p)  =  r 
and  {u  G  ker  DpF{t.  x.p),  Dj:F{t,x,p)u  G  rge  DpF{t,  x,p)}  =>  u  =  0,  respectively,  where 
F{t,x,p)  =  A{t)p  -f  B{t)x  —  b(t),  so  that  they  coincide  with  (i)  and  (ii-a)  of  Theorem  7.1, 
respectively. 
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Remark  7.2.  More  generedly,  it  can  be  shown  that  the  index  of  this  section  always 
coincides  with  the  index  of  [RRh92]  for  the  DAE  (1.1).  The  proof  of  this,  although  concep¬ 
tually  simple,  involves  notational  complications  and  the  detziils  will  not  be  given  here.  But 
we  note  that  the  independence  of  the  index  of  this  paper  from  the  specific  choices  made 
at  each  stage  of  the  reduction  can  be  viewed  as  a  consequence  of  the  fact  that  different 
choices  merely  express  the  geometric,  and  hence  intrinsic,  procedure  of  [RRh92]  in  different 
charts.  □ 

Condition  (ii-b)  is  useful  for  the  study  of  numerical  algorithms:  see  [RRh93]. 

As  a  straightforward  verification  reveals,  conditions  (i)  and  (ii-c)  are  exactly  those  of 
the  case  “jo  =  1”  in  Theorem  3.3  of  [C87]. 

Conditions  (i)  and  (ii-d)  of  Corollary'  7.3  are  used  in  [GrM86]  as  a  criterion  for  “trans¬ 
ferability”,  a  concept  therefore  equivalent  to  the  index  1  condition  of  this  section. 

Acknowledgement:  We  thank  Steve  Campbell  for  pointing  out  to  us  several  contribu¬ 
tions  related  to  this  work  and  for  helpful  discussions. 

References 

[C87]  S.L.  Campbell,  .4  general  form  for  solvable  linear  time  varying  singular  systems  of  differential 
equations,  SIA.M  J  Math.  Anal.  18  (1987),  1101-1115. 

[Co82]  D.  Cobb,  On  the  solutions  of  linear  differential  equations  with  singular  coefficients.  J.  Diff.  Equa¬ 
tions  46  (1982),  311-323. 

[CP83]  S.L.  Campbell  and  L.R.  Petzold,  Canonical  forms  and  solvable  singular  systems  of  differential 
equations,  SIAM  J.  Alg.  Disc.  Meth.  4  (1983),  517-521. 

[D64]  V'.  Dolezal,  The  existence  of  a  continuous  basis  of  certain  linear  subspaces  of  Er  which  depends  on 
a  parameter,  Cas.  Pro.  Pest.  Mat.  89  (1964),  466-468. 

[F65]  K.O.  Friedrich*  Advanced  Ordinary  Differential  Equations.  Gordon  and  Breach,  New  York  (1965). 
[G59]  F.  R.  Gantmacher.  The  Theory  of  Matrices,  Vol  1  and  2,  Translation  from  the  Russian.  Chelsea 
Publ  Co.,  New  York.  NY,  1959. 

[Ge93]  T.  Geerts,  Solvability  for  singular  systems,  Lin.  Alg.  Appl.  181  (1993),  111-130. 

[GrM86]  E.  Griepentrog  and  R.  Maerz,  Differential- Algebraic  Equations  and  Their  Numerical  Treatment. 

vol.  88,  Teubner-Texte  zur  Mathematik,  Leipzig,  1986. 

[H64]  P.  Hartman,  Ordinary  Differential  Equations.  J.  Wiley  and  Sons.  New  York  (1964). 

[K50]  T.  Kato,  On  the  adiabatic  theorem  of  quantum  mechanics.  J.  Phys.  Soc.  Jap.  5  (1950).  435-439. 
[K82]  T.  Kato.  A  Short  Introduction  to  Perturbation  Theory  for  Linear  Operators,  Springer-N’erlag.  New 
York  NY,  1982. 


36 


[KuMe92]  P.  Kunkel  and  \^.  Mehrmann,  Canonical  forms  for  linear  differential-algebraic  equations  with 
variable  coefficients.  Technical  Report  69,  Inst.  f.  Geometrie  und  Prakt.  Mathematik,  RWTH  Aachen 
Germany,  1992. 

[M89a]  R.  Maerz,  Index  2  differential-algebraic  equations.  Results  in  Math.  15  (1989),  149-171 

[.\I89b]  R.  Maerz,  Some  new  results  concerning  index-3  differential-algebraic  equations.  J  Math  Anal 
Appl.  140  (1989),  177-199. 

[RRh91]  P.J.  Rabier  and  W.C.  Rheinboldt,  A  general  existence  and  uniqueness  theory  for  implicit  differ¬ 
ential-algebraic  equations.  Diff.  and  Integral  Equations  4  (1991),  563-582. 

[RRh92]  P.J.  Rabier  and  W.C.  Rheinboldt,  A  geometric  treatment  of  implicit  differential-algebraic  equa¬ 
tions,  J.  Diff.  Equations  (to  appear). 

[RRh93]  P.J.  Rabier  and  W.C.  Rheinboldt,  Finite  differences  for  time  dependent  linear  DAE's.  .Appl 
Math.  Letters  (submitted). 

[S66]  L.  Schwartz,  Theone  des  Distributions.  Hermann,  Paris,  France.  1966. 

[SiB70]  L.M.  Silverman  and  R.S.  Bucy,  Generalizations  of  a  theorem  of  Dolezal.  Math.  Systems  Theory  4 
(1970),  334-339. 

[VLKaSl]  G.C.  Verghese,  B.C.  Levy  and  T.  Kailath,  A  generalized  state-space  for  singular  systems.  IEEE 
Trans.  Automat.  Control  AC-  26  (1981),  811-831. 


37 


REPORT  DOCUMENTATION  PAGE 


.01  i. 'eocft'Og  o^^foen  ":r  :r- s  p'.’.cn  *  •»  '  ■  *.«•' r':*' 

jirenno  <jn<3  'T' %n>t  4ir  inq  {rt|»  a  Jld  niH?0<K3  J^Q  ::iT>oietinq  ir  q ’*»v ••'.»  '  P'"* 

.1  pct  on  .*  Tf^ffT^liOn  *'(  .^0‘nq  ^,1.  '  '  'Pautirq  Ow'OPr  :.  .•,  q*. , -i 

<*.isHqr,^4»  '>«,ie'.’C4  C.''.rqrrp  JPQ  tc  <*'•  _•*  -  :•  V 


1  AGENCY  USE  ONLY  (Leave  oiariK)  2  REPORT  DATE 

’()-  15-q? 


4,  TITLE  AND  SUBTITLE 

CLASSICAL  AND  GENERALISED 
SOLUTIONS  OF  TIME-DEPENDFNT  LINEAR 
DIFFERENTIAL  ALGEBRAIC  EQUATIONS 


6.  AUTHOR(S) 

Pa^'T  i  rk  J  .  Rah  i  er 
Werner  C.  Rheinbnld+ 


7.  PERFORMING  0RGANI2ATI0N  NAME(S)  ANO  AOORESS(ES) 

Deoartmen-^  o'  Mathema+ics  and  S'^a+is+irs 
Universi+v  of  Pi++shurnh 


■CrTi  -cpfOiefJ 

0M8  ‘Vo  0188 


3  REPORT  TYPE  ANO  DATES  COVERED 

tecmnucal  report 


5  FUNDING  NUMBERS 

ONR-N-  onOM-QO-j  -  lOCS 
i^'SF  -CCR  q20T4PA 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  ANO  AOORESS(ES) 

ONR 

NSF 

10.  SPONSORING  MONITORING 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

i 

i 

12a.  DISTRIBUTION  AVAILABILITY  STATEMENT 

Approved  for  public  release  distribution  unlimited 

12b.  DISTRIBUTION  CODE  j 

1 

j 

i 

1 

i 

13.  ABSTRACT  (Maximum  200  woras>  j 

A  reduction  procedure  is  developed  for  linear  + ' me-dependen+  differential  alaehraicj 
equations  (DAFs)  that  transforms  their  solu'''ion5  into  sr^|u+ions  of  smaller  systems  of 
ordi'^ary  differential  equations  (ODFs).  The  procedure  applies  to  '“lassical  as  well 
as  dis+ribijtion  solutions.  In  the  case  of  analytic  coefficients  the  hypotheses 
required  for  the  reduction  no'*'  only  are  necessary  for  the  ''alidi'''y  of  the  existence 
and  uniqueness  results,  hut  they  even  allow  for  the  presence  of  s i nqu  I  ar i t i es . 
Straightforward  extensions  inrludinq  undetermined  systems  and  systems  with  non-ana  I 
non-ana  I yt i c  coefficients  are  also  discussed. 


14.  SUBJECT  TERMS 

linear  taE,  classical  solut  ons,  distribution  solutions 
boundary  value  problems 


IS.  NUMBER  OF  PAGES 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION  20.  LIMITATION  OF  ABSTRACT 
OF  REPORT  OF  THIS  PAGE  OFABSTJ1ACT, 

jjp^l^c^^ifjed  Hnriacstfted  unclassified 


unc I  ass i f i ed 


NSN  7540-01-280-5500 


5taraarc  -orm  298  (Rev  2-39,' 

Pf**Sv'  D>  -NS.  '.Ifl  V39-  H 


GENERAL  INSTRUCTIONS  FOR  COMPLETING  SF  298 


The  Report  Documentation  Page  (RDP)  is  useo  m  announcing  and  cataloging  reports  it  is  important 
that  this  information  be  consistent  with  the  rest  of  the  report,  particularly  the  cover  and  title  page 
Instructions  for  filling  m  each  block  of  the  form  follow  It  is  important  to  stay  within  the  lines  to  meet 
optical  scanning  requirements. 


Block  1  Aqencv  Use  Only  fieave  b/an/c 


Block  2.  Report  Date.  Full  publication  date 
including  day,  month,  and  year,  if  available  (e  g.  1 
Jan  88).  Must  cite  at  least  the  year. 

Blocks.  Type  of  Report  and  Dates  Covered 


State  whether  report  is  interim,  final,  etc.  if 
applicable,  enter  inclusive  report  dates  (e  g.  10 
Jun  87  -  30  Jun  88). 

Block  4.  Title  and  Subtitle.  A  title  is  taken  from 


the  part  of  the  report  that  provides  the  most 
meaningful  and  complete  information.  When  a 
report  is  prepared  in  more  than  one  volume, 
repeat  the  primary  title,  add  volume  number,  and 
include  subtitle  for  the  specific  volume.  On 
classified  documents  enter  the  title  classification 
in  parentheses. 

Blocks.  Fundinq  Numbers.  To  include  contract 


and  grant  numbers;  may  include  program 
element  number(s),  project  number(s),  task 
number(s),  and  work  unit  number(s).  Use  the 
following  labels: 


Contract 

PR  - 

Project 

Grant 

TA  - 

Task 

Program 

WU  - 

Work  Unit 

Element 

Accession  No. 

Blocks.  Author(s).  Name(s)  of  person(s) 
responsible  for  writing  the  report,  performing 
the  research,  or  credited  with  the  content  of  the 
report.  If  editor  or  compiler,  this  should  follow 
the  name(s). 

Block?.  Performinq  Organization  Name(s)  and 


Address(es).  Self-explanatory. 

Block  8.  Performinq  Organization  Report 


Number.  Enter  the  unique  alphanumeric  report 
number(s)  assigned  by  the  organization 
performing  the  report. 

Blocks.  Soonsorinq/Monitorinq  Agency  Name(s) 


and  Address(es).  Self-explanatory. 

Block  10.  Sponsormq/Monitorinq  Agency 
Report  Number  (If  known) 

Block  11.  Supplementary  Notes.  Enter 
information  not  included  elsewhere  such  as: 
Prepared  m  cooperation  with...;  Trans,  of...;  To  be 
published  m....  When  a  report  is  revised,  include 
a  statement  whether  the  new  report  supersedes 
or  supplements  the  older  report. 


Block  12a.  Distribution/Availability  Statement. 


Denotes  public  availability  or  limitations.  Cite  any 
availability  to  the  public.  Enter  additional 
limitations  or  special  markings  m  all  capitals  (e  g. 
NOFORN,  REL,  ITAR), 


See  DoDD  5230.24,  "Distribution 
Statements  on  Technical 
Documents." 

See  authorities. 

See  Handbook  NHB  2200  2. 

Leave  blank 


DOE 

NASA 

NTIS 


Block  12b.  Distribution  Code. 


NASA- 
NTIS  - 


Leave  blank. 

sinter  DOE  distribution  categories 
v'om  the  Standard  Distribution  for 
Unclassified  Scientific  and  Technical 
Reports. 

Leave  blank. 

Leave  blank. 


Block  13.  Abstract.  Include  a  brief  ('Maximum 
200  words)  factual  summary  of  the  most 
significant  information  contained  in  the  report 

Block  14.  Subject  Terms.  Keywords  or  phrases 
identifying  major  subjects  in  the  report. 

Block  15.  Number  of  Pages  Enter  the  total 


number  of  pages 

Block  16.  Price  Code  Enter  appropriate  price 
code  (NTIS  only) 

Blocks  17.  - 19.  Security  Classifications.  Self- 


explanatory  Enter  U.S.  Security  Classification  in 
accordance  with  U.S.  Security  Regulations  (i.e., 
UNCLASSIFIED).  If  form  contains  classified 
information,  stamp  classification  on  the  top  and 
bottom  of  the  page. 

Block  20.  Limitation  of  Abstract  This  block  must 


be  completed  to  assign  a  limitation  to  the 
abstract.  Enter  either  UL  (unlimited)  or  SAR  (same 
as  report).  An  entry  in  this  block  is  necessary  if 
the  abstract  is  to  be  limited.  If  blank,  the  abstract 
IS  assumed  to  be  unlimited. 


StancJard  form  298  Back  (Rev  2-89) 


